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ABSTRACT 

I 

A program  of  research  to  explore  new  display  device  techniques 
and  digital  aystem  architecture  concepts  which  can  be  used  to  realize 
highly  reliable,  low  cost  Interactive  graphic  display  systems  has  been 
completed  at  the  University  of  Illinois  Coordinated  Science  Laboratory. 

* The  techniques  and  concepts  which  were  studied  are  generally  applicable  to 
future  large  scale  computer-based  Information  systems  such  as  those  used 

; in  ground  based  operations  as  well  as  In  specialized  computer-based 

Information  systems  such  as  those  used  In  aircraft,  submarine,  and  shipboard 
Instrumentation  systems.  The  program  during  the  last  contract  period  was 

* divided  into  four  major  areas. 

A major  part  of  the  program  has  been  the  development  of  new 
j techniques  for  the  study  of  the  discharge  processes  In  the  plasma  display 

device,  and  the  development  of  a more  precise  understanding  of  these 
processes.  A second  area  of  Investigation,  based  In  part  on  the  first.  Is 
« the  search  for  new  methods  of  enhancing  the  logical  capability  of  the 

plasma  display  device.  A third  area  Is  an  Investigation  of  spatial  gray 
i scale.  And  finally  a study  has  been  made  of  the  architecture  of  computer- 

based  information  processing  systems  which  emphasize  direct  interaction  of 
^ users  through  graphical  displays. 

- This  report  describes  the  significant  results  in  each  of  the  four 

areas  of  work. 


I 
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I.  INTRODUCTION 


The  diversity  of  topics  discussed  In  this  report  reflects  an 
Interesting  aspect  of  the  present  period  In  plasma  display  development. 

On  the  one  hand  the  availability  of  commercial  devices  with  more  than  250,000 
bistable  addressable  elements  makes  possible  experimental  research  In  the 
architecture  of  computer-based  information  processing  systems,  particularly 
systems  that  emphasize  user  interaction,  and  that  require  graphic  displays. 

On  the  other  hand  the  basic  discharge  processes  on  which  the  operation  of 
the  device  ultimately  depends  are  just  beginning  to  be  understood  In 
detail.  There  Is  thus  both  need  and  opportunity  to  do  fundamental  device 
research.  Both  areas  of  Investigation  are  discussed  in  this  report. 

One  reason  that  more  fundamental  device  research  has  not  been 
done  Is  that  traditional  computer  simulation  of  the  complicated  electrical 
discharge  processes  has  been  expensive.  However,  by  running  computer 
programs  as  part  of  the  background  processing  In  the  PLATO  computer-based 
education  environment  It  has  been  possible  to  simulate  the  discharge 
processes  with  a detail  not  attempted  before.  Another  deterrent  to  basic 
research  has  been  that  experimental  techniques  usually  appropriate  for  gas 
discharge  studies  are  difficult  In  the  environment  of  the  plasma  display 
cell  for  which  a characteristic  dimension  Is  about  30  microns.  New 
techniques  which  Include  electromagnetic  Interactions  from  radio  to  optical 
frequencies  have  been  devised  to  obtain  data  from  the  discharge.  One 
result  of  this  work  Is  the  development  of  a new  interferometric  technique 
for  measuring  phase  shifts  on  smaller  than  10  ^ Angstroms.  This  work  and  its 
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application  to  plasma  cell  measurements  are  described  in  detail  in  the 
following  section  along  with  results  of  both  computer  calculation  and  other 
experimental  studies. 

An  application  of  the  basic  studies  is  the  enhancement  of  logic 
through  control  of  the  gas  discharge  processes.  While  the  preliminary 
results  are  mostly  conceptual,  the  techniques  seem  to  be  sound,  and  the 
importance  of  simplifying  addressing  of  plasma  displays  warrants  discussion 
here. 

The  importance  of  finding  ways  of  encoding  images  for  transmission 
in  information  processing  systems  has  led  to  research  on  algorithms  for 
providing  gray  scale  through  dot  density  control.  In  particular  a simple, 
inexpensive  realization  has  been  found  and  implemented  for  a new  algorithm 
developed  at  Bell  Laboratories.  Studies  on  data  compression  techniques  for 
this  system  are  also  discussed. 

Finally,  systems  studies  have  provided  results  in  two  areas. 

First,  results  from  the  terminal  system/host  system  architecture  studies 
have  provided  definite  guidelines  as  to  how  logical  processing  power  and 
memory  should  be  allocated  in  computer-based  information  and  instrumentation 
systems  which  consist  of  a host  system  and  remote  flat  panel  graphics 
terminals  which  are  connected  via  low  bandwidth  consnvnication  lines.  These 
results  are  applicable  to  both  ground  based  information  systems  as  well  as 
aircraft,  submarine,  and  shipboard  instrumentation  and  information  systems. 

Second,  specific  terminal  architecture  considerations  related  to 
the  results  of  the  first  area  have  emerged.  Microprocessor  based  display 
terminals  and  display  terminals  which  used  considerable  local  processing 
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power  and  memory  located  near  or  on  the  display  have  undergone  extensive 
study.  Specific  recommendations  for  future  terminal  structures  have  been 
made . 

Each  of  the  research  areas  discussed  in  this  report  Is  In  an 
early  phase  of  development.  Each  is  important  to  the  development  of  plasma 
display  technology,  and  we  expect  that  each  will  receive  Increasing 
emphasis  over  the  next  few  years. 


' - t 
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II.  PHYSICS  OF  AC  DISCHARGES 


2.1  Introduction 

This  section  describes  research  that  is  aimed  at  a complete 
understanding  of  the  many  gas  discharge  processes  found  in  the  AC  plasma 
display  panel.  The  program  includes  computer  modeling  and  development  of 
measurement  techniques  needed  to  confirm  the  model.  The  major  results 
show  the  existence  of  a large  plasma  in  the  gas  volume  for  tens  of  micro- 
seconds after  the  discharge.  Neutral  atom  densities  can  be  measured  as  a 
function  of  space  and  time,  and  be  used  to  indicate  the  regions  of  intense 
discharge  activity.  Also  the  experiments  show  that  the  ion  secondary 
effects  mechanism  predominates  over  the  photon  mechanism  for  soda  glass 
cathodes . 

2.2  Computer  Model 

2.2.1  Method  of  Calculation 

Recently,  a number  of  papers  have  appeared  that  model  the  plasma 
panel  discharge  process.  [1,2,3, 4]  All  use  finite  difference  techniques 
with  varying  degrees  of  complexity.  Veron  and  Wang  [1]  used  the  most 
rigorous  method  of  calculation  which  was  originally  developed  by 
Ward  [5, 6, 7, 8, 9].  With  this  very  time  consuming  technique,  the  calculation 
of  the  processes  occurring  for  1 microsecond  in  a plasma  panel  can  take 
about  1 hour  of  CPU  time  on  a large  scale  computer.  Because  of  the  obvious 
cost  factors,  Veron  and  Wang  did  not  run  very  much  data.  Lay  et  al.  [3], 
circumvented  the  CPU  time  problem  by  neglecting  the  field  distortion  in 
the  discharge  gap  caused  by  the  large  densities  of  ions  and  electrons.  As 
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will  be  shown  later  in  this  paper,  this  is  a very  poor  approximation. 
Lanza  [4]  has  developed  a very  clever  technique  for  reducing  the  CPU  time 
required  without  neglecting  field  distortion.  Although  Lanza  presents 
excellent  agreement  between  experimental  and  calculated  currents,  the 
calculations  and  measurements  described  below  tend  to  show  significant 
errors . 

The  basic  calculational  technique  uses  a first  order  finite 
difference  technique  to  solve  the  following  equations  in  one  dimension. 
The  basic  three  equations  to  be  solved  are:  the  electron  continuity 

equation 


Sn  j, 

» 6 3 On  V - -sr—  (n  V ) 
at  ee  ox  e e 


(1) 


which  qualitively  says  the  rate  of  change  of  electron  density  at  a point 
&ne 

in  space  ) i-8  equal  to  the  rate  of  electron  generation  in  the  gap 
(oneV  ) plus  the  rate  of  the  electrons  coming  in  minus  the  rate  of  electrons 
going  out  of  that  point  (-  ^(neve))‘  Similarly  the  ion  continuity  equation 
is 


an  V + -5-  (n.V. ) 
e e Sx  ' i l 


(2) 


where  ne  is  the  electron  number  density,  n^  is  the  ion  number  density, 
is  the  electron  velocity,  is  the  ion  velocity,  a is  the  Townsend 
coefficient  that  describes  the  number  of  electrons  freed  by  the  avalanche 
process  per  unit  length,  x is  space,  and  t is  time.  The  variables  ne>  n^, 
, V^,  and  a are  all  functions  of  both  x and  t.  The  variables  Vg,  V^, 
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and  a are  characterist ic  of  the  gas  and  are  directly  dependent  on  the  value 
of  the  electric  field  E(x,t)  which  must  be  found  from  Poisson's  equation: 


8h 

Sx 


E~  (ni'ne) 


(3) 


where  q is  the  charge  of  the  electron  and  Eq  is  the  dielectric  constant  of 
free  space. 

From  the  above  three  equations  and  knowledge  of  the  functional 
relationship  of  Vgl  V^,  and  cron  E,  one  can  find  the  desired  solution  which 
is  ng  and  n^  for  all  values  of  space  and  time  if  the  necessary  boundary 
conditions  and  initial  conditions  are  specified.  In  a gas  discharge  a 
number  of  physical  effects  can  occur  at  the  cathode  that  will  define  the 
cathode  boundary  conditions.  If  we  define  an  electron  current  density 
Je(x,t)  = n^fx.t)  • Ve(x,t),  then  this  cathode  boundary  condition  is 
specified  by  the  electron  current  from  the  cathode 

Jg(0,t)  - Jo  + ViJ1(0,t)  + fa(x,t)  je(x,t)  dx  (4) 

where  1 is  a very  small  constant  current  emitted  from  the  cathode  to 
get  the  discharge  process  started,  y^  and  y^  are  the  Townsend  coefficients 
for  electrons  emitted  from  the  cathode  due  to  ion  bombardment  and  photon 
bombardment  respectively,  d is  the  distance  between  the  cathode  and  the 
anode  and  a is  a function  of  x and  t that  defines  the  number  of  photons  per 
unit  distance  created  in  an  avalanche  by  an  electron.  The  value  of  o(x,t) 
is  determined  from  E(x,t)  in  much  the  same  manner  as  is  a(x,t).  y . , y 


and  d are  constants  that  characterize  the  discharge. 
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Since  ions  are  not  allowed  to  be  emitted  from  the  cathode 


Jt(0,t)  - 0 


Similarly  nothing  is  allowed  to  be  emitted  by  the  anode,  thus 


Je(d,t)  = 0 


j^d.t)  = 0 


Equations  (4, 5, 6, 7)  define  the  boundary  conditions  that  determine 
ng  and  n^.  To  determine  E(x,t),  we  must  use  Kirkhoff's  voltage  law  and  sum 
around  the  loop  that  includes  the  plasma  display  cell  and  the  applied 
voltage  source  (sustainer).  The  loop  equation  is 


V ‘E  1 + E i + Ed+^  f f(n-n, ) dxdx' 
s c a o €JJ  ei 


where  E and  E are  the  fields  in  the  dielectric  layers  at  the  cathode  and 
the  anode,  l is  the  thickness  of  these  dielectric  layers,  and  Eq  is  the 
field  in  the  gas  volume  at  x=0.  To  find  the  values  of  the  fields  in  the 
dielectric,  E and  E , one  must  know  the  value  of  the  charges  Q and  Q at 

C fl  o J 

the  dielectric  gas  interfaces.  From  the  standard  electric  field  boundary 
conditions  one  gets 


€ E - E - 
r c o € 
o 


E . - € E - 
d r a € 

o 


and 
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where  is  the  relative  dielectric  constant  of  the  dielectric  layers  and 

Ej  is  the  electric  field  in  the  gas  volume  at  x=d. 

The  charges  Qq  and  are  obtained  by  integrating  all  currents 

flowing  into  and  out  of  the  cathode  and  anode. 

For  initial  conditions  one  must  specify  ne(x,0),  n^(x,0),  QQ(t”0), 

and  Q{j(t=0).  The  sustain  voltage  Vg  is  the  desired  function  of  time. 

The  gas  parameters  used  those  given  by  Ward  [5].  The  dependence 

of  a,  V , and  V.  on  the  E field  is 
e i 

a *=  pA  expt-BCp/E)1^2]  (11) 


where  p is  the  gas  pressure  and  A and  B are  constants  characteristic  of 
the  gas. 


p E 

V - — — 
e p 


(12) 


where  p,  is  the  electron  mobility. 


U+E 


CE, 


i p E/15W1 

k E1/2 

„ „ + /,  ric,-3/2  3/2, 

\/2  ^ P 


E/p  > w. 


(13a) 

(13b) 


where  u+,  k+,  C,  D,  and  w^  are  constants  characteristic  of  the  gas.  Ward  [5] 
gives  values  of  these  constants  for  various  gases.  In  the  examples  used 
here,  pure  Ne  was  used  and  has  the  following  constants 

A ■ 8.2  cm  * torr  1 

_ „ ,A/2  -1/2  „ -1/2 

B"  17.0  volt  cm  torr 


-3  -1 

C « 6 x 10  torr  cm  volt 
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n = / 1 q q ,3/2  -3/2  _ -3/2 

D = 41.98  volt  cm  torr 

2 -1-1 
li  = 10  cm  torr  volt  sec 


32  -1  -i 

= 2.9  x 10  cm  torr  volt  sec 

. . , ,_4  3/2  „ 1/2  , -1/2  -1 

k+  * 1.7  x 10  cm  torr  volt  sec 


The  calculational  procedure  follows.  In  finite  difference  calcu- 
lations, each  iteration  produces  a solution  advanced  by  a time  increment  At. 
This  value  is  chosen  from  a convergence  criteria  [9]: 

^ <max|ve(x,t)|x  (14) 

where  Ax  is  the  distance  between  mesh  points.  This  relation  simply  states 
that  one  must  calculate  the  next  iteration  before  the  fastest  electrons 
have  time  to  get  from  one  mesh  point  to  the  next.  Failure  to  obey  (14) 
causes  the  calculation  to  blow  up.  For  30  mesh  points  across  the  discharge 
gap  and  for  typical  plasma  panel  conditions,  At  is  about  10  ps . Thus  for 
1 (js  of  plasma  panel  data,  the  following  sequence  is  iterated  10"*  times. 


1.  Given  the  initial  valuer,  of  Qq  , Q^,  V^(t),  n^(x,t),  n^(x,t),  calculate 
for  all  x E(x,t)  from  (3),  (8),  (9),  and  (10). 

2.  Calculate  Vg(x,t)  * — E(x,t)  for  all  x. 

3.  Calculate  for  all  x j (x,t)  * V (x,t)  n (x,t). 

4.  Calculate  V^(x,t)  from  (13)  for  all  x. 

5.  Calculate  J^(x,t)  ■ V^x.t)  n^(x,t)  for  all  x. 

6.  Determine  At  ■ ^ i ■ — 7-1—  . 

max|Ve(x,t)|x 

7.  Calculate  t ■ t + At. 

Calculate  ne(x,t)  from  a difference  equation  derived  from  (1)  with  the 
aid  of  (4)  for  all  x. 


* 


8. 
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9.  Calculate  n^(x,t)  from  a difference  equation  derived  from  (2). 

10.  Calculate  Q = Q + At  • j (0,t) 

o o i 

V4>e) 

11.  Return  to  step  1. 

2.2.2  Results  of  Calculations 

The  data  presented  here  are  a typical  sample  from  more  than  1000 
hours  of  CPU  computations  performed  over  the  past  2 years.  These  data  were 
computed  in  a background  mode  on  the  CDC  Cyber  73  (6500)  used  by  the  PLATO 
terminals  having  a 512  x 512  AC  plasma  panel.  The  first  3 figures  presented 
here  are  taken  directly  from  a PLATO  terminal.  The  ability  to  rapidly 
observe  and  manipulate  the  various  desired  functions  gives  the  researcher 
a very  useful  qualitative  understanding  of  the  processes  that  has  led  to 
a number  of  creative  advances. 

Figure  2.1  shows  the  applied  voltage  and  discharge  current  pulse 
for  a typical  one  dimensional  plasma  display  cell.  This  is  qualitatively 
similar  to  all  previously  presented  data.  Figure  2.2  shows  the  electron 
charge  density  as  a function  of  the  space  between  the  cathode  and  the  anode. 
This  shows  the  early  exponential  growth  on  the  lower  right  hand  side,  the 
build  up  to  a very  large  electron  density  that  fills  the  entire  gap  except 
for  the  cathode  region,  and  then  the  slow  decay  of  the  large  electron 
density  toward  the  anode.  Each  complete  line  shown  in  Figure  2.2  is  composed 
of  14  straight  line  segments  that  connect  15  data  points  for  15  positions 
in  the  gap.  Thus  the  decay  of  electrons  is  not  accurately  represented  in 
Figure  2.2  due  to  undersampling.  The  peak  of  the  discharge  current  shown 


Figure  2.2.  Log  plot  of  electron  charge  density  as  a 
function  of  space  and  time.  The  cathode 
is  on  left  and  the  anode  on  right.  A time 
increment  of  10  ns  separates  each  curve. 
The  origin  of  the  vertical  axis  is  10"^® 
coulombs /cm^ . 


13 


In  Figure  2.1  closely  corresponds  to  the  curve  showing  the  maximum  total 
electron  density  in  Figure  2.2. 

The  major  difference  between  previous  calculations  and  those 
presented  here  is  the  very  large  and  slowly  decaying  electron  density 
shown  in  Figure  2.2.  This  is  typical  of  all  data  computed  here  using 
Ward's  finite  difference  technique.  If  one  looks  at  the  corresponding 
ion  density  distribution,  one  finds  it  almost  identical  to  the  electron 
density  after  the  peak  of  the  discharge  current.  This  is  in  fact  a plasma 
as  defined  by  the  strict  sense  that  the  Debye  length  is  much  shorter  than 
the  cell  dimensions.  The  electron  and  ion  densities  are  so  great  that 
they  cancel  out  most  of  the  electric  field  on  the  anode  side  of  the  gap. 
With  little  electric  field,  the  electrons  and  ions  can  move  only  by  the 
slow  diffusion  process  or  can  slowly  recombine  and  are  thus  effectively 
trapped.  Typically  this  plasma  region  exists  for  a number  of  microseconds 
after  the  discharge.  It  is  not  uncommon  for  the  plasma  to  still  exist 
when  the  sustain  waveform  falls  5 or  10  microseconds  after  the  discharge. 
Thus  with  certain  sustain  waveforms,  the  existence  of  the  plasma  from  the 
previous  discharge  can  be  a very  important  initial  boundary  condition. 

The  electric  field  as  a function  of  gap  spacing  is  shown  in 
Figure  2.3.  The  field  is  initially  equal  for  every  point  along  the  gap 
and  increases  as  the  sustain  waveform  rises.  When  the  discharge  occurs, 
the  field  at  the  anode  decreases  due  to  the  growth  of  the  plasma,  and 
thus  makes  the  field  at  the  cathode  very  large.  It  is  not  uncommon  to 
see  cathode  fields  larger  than  10^  volts/cm.  It  is  interesting  to  note 
that  field  emission  has  been  observed  for  these  fields  intensities  and 
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thus  might  be  an  important  cathode  emission  process  at  the  peak  of  the 
discharge.  Note  also  that  most  of  the  discharge  activity  or  light  output 
occurs  near  the  cathode  which  is  consistent  with  Anderson's  measurements  [10]. 
Also  the  field  in  the  plasma  region  is  near  zero. 

2.3  Capacitance  Measurements 

The  existence  of  the  plasma  in  the  computer  model  stimulated  the 
realization  that  the  capacitance  of  the  plasma  cell  should  be  varying  as 

a function  of  the  plasma.  Figure  2.2  shows  the  density  of  the  plasma  is 

"63  13  3 

about  10  C/cm  or  about  10  electrons/cm  . This  corresponds  to  a 

plasma  frequency  of  about  50  GHz.  Thus  for  frequencies  considerably  below 

this,  the  plasma  will  appear  as  a material  having  a near  infinite  dielectric 

constant.  This  will  change  the  capacitance  of  the  cell  as  shown  in  Figure 

2.4.  Before  the  discharge,  there  is  no  plasma  and  the  capacitance  is  the 

neutral  value  CQ.  At  the  discharge  time  t^,  the  cell  is  filled  with  plasma 

and  the  capacitance  is  at  a maximum  value  C^.  As  the  plasma  decays  the 

capacitance  decreases  until  the  plasma  is  gone  and  the  initial  Cq  value  is 

restored.  Thus  a measurement  of  the  capacitance  of  a cell  gives  information 

about  the  spatial  extent  and  the  existence  of  the  plasma.  Note  that  at 

the  frequencies  where  such  capacitance  measurements  are  conveniently  made, 

no  information  is  obtained  about  the  plasma  density  since  these  frequencies 

are  much  lower  than  the  plasma  frequency. 

Figure  2.5  shows  a simple  method  of  measuring  a capacitance 
change  in  a plasma  cell.  It  is  necessary  that  the  capacitance  measuring 
circuit  elements  not  disturb  the  application  of  the  sustain  voltage  to  the 
cell.  Also  the  sustain  voltage  signal  should  not  cause  excess  noise  to 
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obscure  the  capacitance  measurement.  For  these  reasons  It  Is  usually 
convenient  to  use  a rf  signal  In  the  range  of  10  to  100  MHz.  In  the  Figure 
2.5  circuit,  a resonance  Is  achieved  between  the  variable  Inductor  and  the 
capacitance  of  the  plasma  cell.  The  capacitor  is  large  compared  to  the 
cell  capacitance  but  small  enough  so  as  not  to  overload  the  sustain 
generator.  The  variable  inductor  must  be  small  enough  so  as  not  to  cause 
significant  ringing  when  the  sustain  generator  rises.  The  oscillator  and 
the  oscilloscope  are  loosely  coupled  to  the  tuned  circuit  by  the  small 
inductors.  This  technique  is  very  sensitive  to  small  capacitance  changes 
if  the  Q of  the  LC  circuit  is  very  high.  Thus  it  is  particularly  useful 
for  measuring  the  capacitance  changes  in  standard  plasma  panels  where  the 
small  cell  capacitance  is  in  parallel  with  the  many  stray  electrode 
capacitances  making  the  relative  capacitance  change  very  small. 

Figures  2.6  and  2.7  show  the  discharge  current  and  the  rf 
envelope  measured  with  the  circuit  of  Figure  2.5  on  a standard  128-33 
Digivue  panel.  The  envelope  waveform  is  strongly  dependent  on  the 
frequency  of  the  oscillator.  Figure  2.6  shows  the  oscillator  tuned  to 
resonate  with  the  neutral  capacitance  value  Cq.  When  the  discharge  occurs, 
the  increased  capacitance  tunes  the  circuit  out  of  resonance  with  the 
oscillator.  At  a slightly  different  oscillator  frequency,  Figure  2.7 
shows  the  changing  capacitance  tuning  the  circuit  in  and  out  of  resonance. 
Note  that  for  this  particular  standard  plasma  panel,  the  capacitance  is 
still  changing,  15  microseconds  after  the  discharge.  This  is  excellent 
experimental  proof  of  the  existence  of  the  plasma  predicted  by  the  computer 
calculations  shown  in  Figure  2.2. 
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Figure  2.6.  RF  envelope  of  19  Hllz  signal  in  upper  trace. 

Discharge  current  at  40  ma  per  division  in 
lower  trace.  Horizontal  scale:  2 p.s  per 

division.  Data  taken  with  circuit  of  Figure 
2.5  on  standard  128-33  Digivue  panel. 
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The  circuit  of  Figure  2.5  has  a few  disadvantages.  It  is  noisy 

since  the  rise  of  the  sustainer  always  causes  a considerable  amount  of 

\ 

ringing.  This  may  die  out  in  a few  microseconds,  but  it  obscures  the 
discharge  region.  A second  disadvantage  is  the  complex  relationship 
between  the  envelope  amplitude  and  the  capacitance  that  depends  strongly 
on  circuit  Q and  oscillator  frequency.  To  solve  these  problems,  the  circuit 
of  Figure  2.8  was  used.  Quarter  wave  stubs  are  used  to  isolate  the 
sustainer  from  the  plasma  cell  at  the  rf  frequency.  The  capacitors  of 
value  C^,  are  large  enough  to  act  as  a short  circuit  at  the  rf  frequency, 
but  not  overload  the  sustain  generator.  The  final  quarter  wave  stub,  at 
the  oscilloscope,  effectively  filters  out  the  sustain  generator  frequencies 
but  lets  the  rf  frequency  pass.  Thus  at  the  rf  frequency  the  circuit 
behaves  as  if  one  lead  of  the  cell  is  connected  to  the  rf  generator  and 
the  other  to  the  oscilloscope.  Because  of  oscilloscope  capacitance,  the 
measured  signal  is  a voltage  division  between  this  capacitance  and  the 
cell  capacitance.  With  the  knowledge  of  the  oscilloscope  probe  capacitance, 
the  actual  cell  capacitance  can  be  easily  calculated. 

To  reduce  the  problem  of  the  stray  capacitances  that  usually 
shunt  out  the  desired  cell  capacitance,  a very  large  area  cell  having 
disc  shaped  electrodes  1 cm  in  diameter  was  constructed.  The  dielectric 
glass  layer  was  about  25  pm  thick  and  the  cathode -anode  gap  spacing  a bout 
200  pm.  Figure  2.9  shows  the  output  waveform  obtained  from  this  special 
cell  by  means  of  the  circuit  of  Figure  2.8.  Note  that  the  plasma  exists 
for  almost  60  pa  after  the  discharge.  This  time  varies  with  gas  mixture 
and  pressure.  The  longest  time  observed  was  100  ps.  As  a general  rule 
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Figure  2.8.  Circuit  diagram  of  quarter  wave  stub  technique 
for  measuring  the  cell  capacitance  variation. 


RF  envelope  of  57  MHz  signal  in  upper  trace. 
Sustain  voltage  waveform  at  50  volts  per 
division  in  lower  trace.  Horizontal  scale: 

20  p,s  per  division.  Data  taken  with  circuit 
of  Figure  2.8  on  special  1 cm  diameter  plasma 
cell  containing  250  torr  Ne  and  0.25  torr  Xe. 
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neon-xenon  mixtures  have  a considerably  longer  plasma  lifetime  than 
neon-argon  mixtures.  This  can  be  attributed  to  the  diffusion  and  recombi- 
nation coefficients  of  the  minority  gas  ions. 

The  actual  capacitance  values  of  the  cell  can  be  determined  from 
the  data  shown  in  Figure  2.9.  This  again  varies  with  gas  pressure  and 
mixture,  but  typical  values,  referring  to  Figure  2.4,  are  7 pf  for  Cq  and 
50  pf  for  0^.  From  these  numbers  one  can  estimate  that  at  the  peak,  the 
plasma  fills  more  than  6/7  of  the  cathode-anode  gap.  This  is  in  agreement 
with  calculations  such  as  those  shown  in  Figure  2.2. 

2.4  Microwave  Measurements 

Although  the  capacitance  is  a relatively  easy  quantity  to 
measure,  it  gives  virtually  no  information  about  the  density  of  the  plasma. 
For  this  one  must  go  to  the  higher  frequencies.  This  means  microwave  or 
optical  frequencies.  Because  the  wavelength  of  the  microwaves  is  very 
much  larger  than  typical  plasma  cell  dimensions,  a very  special  cell 
structure  must  be  used.  Figure  2.10  shows  a plasma  cell  that  acts  as  a 
microwave  strip  transmission  line.  The  impedance  of  this  line  is  about 
50  ohms.  The  sustain  generator  is  connected  between  the  ground  plane 
electrode  and  transmission  electrode  by  means  of  a very  thin  high  impedance 
wire.  10  GHz  microwaves  were  sent  down  the  strip  line,  and  the  output  was 
detected  to  show  the  transmission  characteristics.  Figure  2.11  shows 
typical  microwave  transmission  data  along  with  simultaneously  measured 
low  frequency  capacitance  measurements.  The  microwave  transmission  is 
very  effectivel'  cut  off  after  the  discharge.  After  20  us  the  transmission 
begins  to  pickup  until  it  has  nearly  perfect  transmission  after  50  us- 
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Fgiure  2.10.  Microwave  stripline  plasma  cell. 
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Figure  2.11.  Microwave  transmission  data  in  lower  trace  and 
low  frequency  capacitance  envelope  in  upper 
trace,  for  strip  line  cell  shown  in  Figure  2.10, 
containing  250  torr  Neon  and  .4  torr  Argon. 
Horizontal  scale  is  20  p,sec  per  division. 
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The  interpretation  of  this  data  is  complicated  since  the 

transmission  coefficient  should  be  dependent  on  both  the  plasma  density 

and  its  spatial  extent  throughout  the  gap.  However  one  might  roughly 

estimate  that  when  the  transmission  is  cut  off  the  plasma  density  is 

great  enough  so  that  its  plasma  frequency  is  greater  than  the  10  GHz 

microwave  frequency.  This  would  indicate  that  the  density  is  greater 
12  '3 

than  10  cm  . However  this  estimate  is  further  complicated  by  the  fact 
that  absorptive  losses  due  to  the  plasma  collisions  with  the  neutral  gas 
atoms  can  be  very  significant  at  the  pressures  typically  used  for  plasma 
panels.  Thus  much  additional  data  at  other  microwave  frequencies  is 
needed . 

2.5  Optical  Interferometry 

One  attractive  method  of  measuring  electron  density  as  a function 
of  time  and  space  in  the  AC  plasma  display  panel  is  optical  interferometry. 
The  basic  principal  is  that  the  presence  of  electrons  in  a volume  will 
when  solving  Maxwell's  equations,  cause  electromagnetic  radiation  to 
propagate  at  a different  velocity.  For  optical  frequencies  this  is 
traditionally  recognized  as  an  index  of  refraction  change.  When  light  is 
sent  through  the  plasma  all  of  a given  length,  the  light  will  experience 
a slight  phase  perturbation  that  can  be  measi'~ed  by  an  interferometer. 

The  optical  Interferometry  method  of  measuring  electron  density 
has  a number  of  advantages  over  other  techniques.  It  usually  does  not 
disturb  the  plasma  such  as  do  electric  probes  or  microwaves.  It  can  make 
very  fast  temporal  measurements  usually  only  limited  by  optical  detector 
speed.  Optical  beams  can  be  focused  to  a small  point  whose  diameter  is 
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limited  by  diffraction  theory  and  thus  they  offer  excellent  spatial 
resolution.  Microwaves,  on  the  other  hand,  have  a wavelength  very  large 
compared  to  a plasma  cell  and  thus  offer  practically  no  spatial  resolution 
at  all. 


2.5.1  An  Accurate  Interferometer 

The  prime  question  that  needs  answering  is  does  the  electron 
density  in  a plasma  panel  have  sufficient  amplitude  to  be  measurable  by 
interferometry  techniques?  From  Maxwell's  equations  we  get  the  well 
known  relations  that 


r 


1 


(15) 


where  is  the  relative  dielectric  constant,  n is  the  index  of  refraction, 
uip  is  the  electron  plasma  frequency,  and  ui  is  the  frequency  at  which  the 
measurement  of  €r  is  made,  in  this  case  the  frequency  of  the  probing  light 
beam.  It  is  also  well  known  that  u>p  varies  as  the  electron  density 
such  that 


uup 


(16) 


where  n^  is  the  electron  density,  €q  is  the  dielectric  constant  of  free 

space,  e is  the  electron  charge  and  m^  is  the  electron  mass.  For  optical 

frequencies  and  the  densities  encountered  in  this  work  u>  » cup  so  that 

the  following  approximation  can  be  made 

. 2 
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from  (16)  and  (17)  it  is  seen  that  for  a given  u),  An  is  proportional  to 
Ane>  Now  if  the  light  beam  traverses  a path  of  length  L,  the  phase  shift 
expressed  in  radians  is 


Acp  - 2IlL(An  - l)/\  (18) 

where  X is  the  wavelength  of  light.  From  (16),  (17),  and  (18)  we  find 
that 

Acpe  - - BAneL\  (19) 

e2  -15 

where  B “ r “ 2.817  x 10  meters,  c being  the  speed  of  light  in 

4H€  C 
o 

vacuum. 

Now  from  the  numerical  calculations  of  Section  2,  we  found 
that  the  peak  electron  density  should  be  in  the  range  of 

12  13  -3 

10  < n < 10  cm  J 

e 

If  we  assume  a very  long  plasma  display  cell  of  1 cm  and  a typical  laser 

o 12  -3 

wavelength  of  6328  A,  from  (19)  and  an  electron  density  of  10  cm  we 

get  a phase  shift  of  1.78  x 10  5 radians  or  0.018  angstroms.  This  is  a 

very  small  phase  change  that  looks  very  discouraging  since  it  is  generally 

considered  that  good  interferometry  technique  can  only  measure  fringe 

shifts  of  Yoo»  [ll].  The  basic  problem  that  limits  the  sensitivity  of 

most  interferometers  is  mechanical  stability.  The  mechanical  elements 

usually  move  many  angstroms  due  to  thermal  expansion  and  sonic  vibrations. 

There  are  two  factors  that  created  great  optimism  in  being 

able  to  measure  the  required  small  fringe  shifts.  The  first  was  the 

use  of  a particular  type  of  interferometer  in  which  most  mechanical 
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motion  was  seen  by  both  legs  of  the  Interferometer,  thus  canceling  out 
to  a relative  motion  of  zero.  The  second  factor  was  that  since  the 
plasma  display  panel  operates  at  high  frequencies  (50  KHz)  relative  to 
the  acoustic  vibrations,  much  of  this  acoustic  noise  can  be  filtered  and 
also  averaging  techniques  are  convenient  to  use  thus  greatly  increasing 
the  signal  to  noise  ratio. 

The  Jamin  interferometer  arrangement  used  in  this  work  is  shown 
in  Figure  2.12.  The  two  optical  flats  are  aluminized  on  their  back  sides 
to  form  mirrors.  The  beam  from  the  laser  is  initially  split  by  the  surface 
reflections  of  the  glass.  The  detector  used  in  these  experiments  was  a 
high  speed  photodiode.  This  interferometer  is  very  simple  in  that  there 
are  only  two  elements  that  can  move  to  cause  an  unwanted  phase  change  - 
the  two  flats.  Translational  motion  of  the  flats  in  the  y and  z directions 
causes  no  change  in  optical  path  length  for  perfect  optical  flats.  For 
motion  in  the  x direction,  both  the  light  beam  going  through  the  plasma 
cell  and  the  reference  beam  experience  the  same  path  change  and  thus  their 
relative  phases  remain  unchanged.  The  only  motion  that  can  cause  a noise 
signal  is  a rotational  motion  about  the  z axis.  Thus  this  is  an  extremely 
stable  interferometer.  Most  other  interferometers  such  as  the  Michelson 
require  a very  stable  mechanical  mount  to  get  the  output  fringes  to 
remain  stationary  enough  for  the  eye  to  recognize  them.  The  Jamin 
interferometer  meets  this  criteria  when  the  flats  are  held  with  the  human 
hand . 

With  the  mechanical  noise  problem  essentially  solved,  other 
noise  sources  become  important.  The  light  source  becomes  the  main 


problem.  Lasers  are  generally  very  monochromatic  with  bright  and  well 
collimated  beams  so  that  they  are  the  only  possible  choice  of  light 
source  for  this  problem.  However,  they  are  also  fairly  noisy  devices.  A 
17.  noise  level  may  seem  good  to  a laser  manufacturer,  but  for  this  problem 
it  represents  a signal  equivalent  to  a \/200  phase  shift  and  thus  is  quite 
unacceptable.  Fortunately  this  noise  can  be  canceled  out  by  using  a second 
detector  and  a differential  amplifier.  The  second  detector  can  be  placed 
at  a number  of  places  along  the  beam.  For  instance  in  Figure  2.12  a beam 
splitter  can  be  placed  between  the  laser  and  the  interferometer.  The 
deflected  beam  from  this  beamsplitter  goes  to  the  second  detector.  The 
outputs  of  both  detectors  are  fed  to  a differential  amplifier.  If  the 
two  detectors  are  sufficiently  linear  and  properly  frequency  compensated, 
all  of  the  excess  laser  noise  can  be  canceled  out  to  the  level  of  the  shot 
noise  plus  the  detector  noise.  The  shot  noise  is  not  the  shot  noise 
traditionally  associated  with  the  leakage  current  of  the  detector,  but 
rather  the  shot  noise  due  to  the  current  from  the  laser  light.  This  noise 
current  is  given  by  the  familiar  formula 

I - \/2 Af  Ie  (20) 

sn 

where  Af  is  the  bandwidth,  I is  the  signal  current  and  e is  the  electron 
charge.  Since  the  signal  we  wish  to  measure  clearly  increases  propor- 
tionally to  the  laser  light  intensity,  the  signal  to  shot  noise  ratio 
increases  as  the  square  root  of  the  laser  intensity.  Thus  it  pays  to  use 
the  largest  power  laser  available.  A high  power  laser  also  eases  the 
electrical  noise  problem  since  the  electrical  noise  such  as  that  due  to 
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the  high  voltage  sustain  waveform  is  not  dependent  on  laser  power.  Thus 
the  signal  to  electrical  noise  ratio  is  proportional  to  the  laser  power. 

With  the  power  from  even  a 1 milliwatt  laser,  the  noise  generated  in  the 
detector  such  as  resistive  noise  is  so  small  compared  to  the  signal  shot 
noise  that  all  other  noise  sources  can  be  neglected. 

There  are  two  basic  improvements  to  the  system  of  Figure  2.12 
that  improve  the  signal  to  noise  ratio  and  the  stability.  To  increase 
stability,  a glass  place  controlled  by  a servo  system  is  placed  in  one  leg 
of  the  interferometer.  A change  in  the  angle  of  the  glass  plate  relative 
to  the  laser  beam  will  cause  the  beam  to  travel  a different  length  path 
through  the  glass.  This  is  a very  old  trick  that  allows  one  to  change 
the  relative  phase  of  one  beam  over  the  other  and  thus  position  the 
output  fringes  at  the  proper  spot.  If  the  flats  are  parallel  to  each  other, 
the  aperture  of  the  laser  beam  will  allow  only  part  of  one  fringe  to  show 
up  at  the  detector  making  the  detector  output  approximate  the  amplitude 
of  a single  point  in  space  while  the  fringes  move  by.  This  amplitude  shows 


up  as  the  familiar  cos*  (II  t)  relation  where  A is  the  optical  path  difference 

between  the  two  interferometer  beams.  Since  the  fringe  movement  measured 

in  these  experiments  is  very  much  smaller  than  \,  the  relationship  between 

the  phase  change  and  the  detector  output  can  be  assumed  linear.  The 

2 

proportionality  constant  is  determined  by  where  along  the  cos  curve  one 
operates.  The  angular  position  of  the  glass  plate  determines  this  operating 

point.  For  maximum  signal  to  noise  ratio,  this  operating  point  should  be 

2 2 A 

where  the  cos  curve  has  maximum  slope.  This  occurs  whenever  cos  (II  — ) * 1/2. 

h 

This  is  also  the  operating  point  where  the  linearity  approximation  holds 
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the  best.  It  should  be  noted  that  the  sign  of  the  proportionality  constant 
can  be  either  positive  or  negative,  depending  on  which  preferred  operating 
point  is  chosen.  This  fact  serves  as  a helpful  check  for  determining  if 
the  measured  signal  is  a real  phase  change  or  some  sort  of  electrical 
noise.  Since  the  phase  change  signal  can  be  inverted  simply  by  slightly 
moving  the  glass  plate,  and  since  any  electrical  noise  is  most  probably 
independent  of  the  position  of  the  plate,  the  noise  and  the  signal  can 
very  easily  be  separated. 

As  discussed  above,  the  Jamin  interferometer  is  very  stable; 
however  over  a period  of  minutes,  small  temperature  changes  can  cause  it 
to  slightly  drift.  Also  room  air  currents  passing  through  the  inter- 
ferometer can  cause  one  leg  to  have  a slightly  different  optical  length 
than  the  other.  Most  of  these  effects  have  frequency  components  that 

are  much  lower  than  those  of  interest  for  this  work,  however,  these 

2 

disturbances  change  the  operating  point  along  the  cos  curve  and  thus 
modulate  the  signals  at  the  higher  frequencies.  For  instance,  if  signal 
integration  techniques  are  used  to  improve  the  signal  to  noise  ratio,  a 
number  of  minutes  might  be  needed  for  the  integration.  Any  drift  of  the 
operating  point  due  to  thermal  change  would  cause  a false  signal  to  appear 
after  the  integration.  Thus  a servo  system  is  needed  to  hole  the  inter- 
ferometer to  a fixed  operating  point.  As  shown  in  Figure  2.13,  the  glass 
plate  is  connected  to  a deflection  motor  which  in  this  case  was  a simple 
speaker  coil  arrangement.  The  coil  was  driven  by  a servo  amplifier  with 
one  differential  input  connected  to  the  detector  output  and  the  other  to 
a variable  voltage  source.  When  operating  properly,  the  glass  plate 
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rotates  until  the  voltage  output  from  the  detector  is  equal  to  the  variable 
voltage  source.  Thus  by  adjusting  this  voltage  source,  one  can  determine 
the  operating  point  of  the  interferometer. 

The  deflection  motor-glass  plate  arrangement  had  a natural 
resonance  at  120  Hz.  The  servo  amp  had  a computer  aided  designed  5 pole 
RC  low  pass  filter  with  a rolloff  very  close  to  12  Db  per  octave.  Since 
oscillations  in  the  servo  system  are  undesirable,  the  gain  of  the  amplifier 
was  adjusted  so  that  the  gain  at  the  120  Hz  resonance  was  just  under  1. 

Thus  the  frequency  response  of  the  servo  system  does  not  affect  the  desired 
signals  which  all  have  frequency  components  above  1 KHz,  but  is  fast  enough 
to  cancel  out  variations  due  to  rapidly  moving  air  currents. 

A second  improvement  over  Figure  2.12  is  the  addition  of  beam- 
splitters on  the  front  surface  of  the  flats  as  shewn  in  Figure  2.13.  This 
serves  to  increase  the  signal  to  noise  ratio.  From  the  above  discussions, 
it  is  clear  that  when  shot  noise  is  limited,  the  signal  to  noise  ratio  will 
be  improved  the  greater  the  signal  at  the  detector.  Also  the  signal  at 
the  detector  should  be  fully  modulated  by  the  interferometer.  Another  way 
of  saying  this  is  that  when  the  interferometer  is  adjusted  for  a minimum, 
the  input  to  the  detector  should  be  zero.  In  Figure  2.13  the  laser  beam  is 
initially  split  by  the  reflectivity  of  glass  surface.  For  the  proper 
polarization  and  a flat  at  45°  to  the  beam,  the  reflectivity  of  the  glass 
surface  is  approximately  107,.  Thus  the  two  beams  are  not  of  equal  intensity. 
When  the  two  beams  are  recombined,  the  stronger  beam  that  goes  through  the 
plasma  cell  only  experiences  a 107,  reflection,  and  thus  most  of  the  laser 
energy  does  not  reach  the  detector.  What  does  reach  the  detector  is  1007 
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modulated,  but  Its  maximum  intensity  can  only  be  about  207.  of  the  laser  beam 
intensity.  The  addition  of  307.  reflecting  beamsplitters  at  the  points  of 
surface  reflection  as  shown  by  the  darkened  areas  on  the  front  of  the  flats 
in  Figure  2.13  neatly  solves  this  problem.  This  makes  both  beams  of  nearly 
equal  intensity.  The  output  to  the  detector  is  1007.  modulated  and  the 
maximum  signal  to  the  detector  can  be  nearly  1007.  of  the  laser  signal. 
Hovever,  as  stated  above  the  interferometer  is  properly  biased  when  the 
ou^'il  is  at  half  maximum  so  that  only  half  of  the  laser  energy  goes  to 
detector  1.  It  is  convenient  to  place  a second  detector  where  the 
remainder  of  the  energy  goes.  This  is  designated  detector  2.  If  the 
outputs  of  both  of  these  detectors  are  fed  into  a differential  amplifier, 
the  laser  noise  canceling  function  can  be  performed.  Detector  2 essentially 
does  the  job  that  a beamsplitter  and  detector  arrangement  between  the 
laser  and  the  interferometer  did  as  discussed  above.  Such  a beamsplitter 
would  have  to  deflect  energy  out  of  the  interferometer  and  thus  reduce  the 
signal  to  noise  ratio.  The  arrangement  of  Figure  2.13  allows  the  total 
laser  power  to  enter  the  interferometer,  along  with  the  added  advantage 
that  the  laser  noise  canceling  detector  also  has  a signal.  With  proper 
care,  the  signals  entering  the  two  detectors  will  be  180°  out  of  phase, 
whereas  the  laser  noise  will  be  in  phase,  thus  the  output  of  the  differential 
amplifier  will  show  only  the  signal  plus  the  shot  noise.  The  balancing  of 
the  signals  to  the  detectors  is  achieved  by  adjusting  the  glass  plate 
position.  For  ideal  beamsplitters  with  507,  reflection  and  transmission, 
tin  two  output  beams  will  be  exactly  balanced  at  the  correct  bias  point. 
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The  lens  shown  in  Figure  2.13  is  used  to  focus  the  laser  beam  in 
the  area  of  the  plasma  cell  to  achieve  optimum  spatial  resolution.  The 
length  of  the  plasma  cell  was  chosen  to  be  about  1 cm,  which  allows  a 
reasonable  compromise  between  the  signal  size  and  the  desired  spatial 
resolution.  The  longer  the  cell,  the  larger  the  signal  according  to 
equation  (19).  However,  diffraction  theory  limits  the  length  over  which 
the  beam  can  stay  focused.  For  the  beam  to  stay  focused  over  1 cm,  the 
beam  waist  had  to  be  made  about  25  nm  in  diameter.  A shorter  length 
would  allow  a smaller  beam  waist.  However  the  25  pm  resolution  seems 
initially  satisfactory  for  the  150-200  um  gap  spacings  used  in  typical 
plasma  cells. 

The  system  shown  in  Figure  2.13  was  developed  specifically  for 
measuring  these  small  phase  shifts.  The  two  optical  flats  define  the  basic 
Jamin  interferometer  which  has  the  desirable  property  that  translational 
motion  of  the  flats  changes  the  length  of  both  optical  paths  equally,  thus 
making  the  output  signal  invariant  to  mechanical  vibrations.  Each  of  the 
flats  has  a small  dielectric  beamsplitter  on  the  front  surface  so  that  the 
two  optical  paths  have  equal  intensity.  Thus  a half  wavelength  phase 
change  in  one  path  can  cause  100?!.  modulation  of  the  output  beam.  A glass 
plate  is  inserted  in  one  optical  path  and  is  rotated  by  a servo  system  to 
make  slight  corrections  for  low  frequency  thermal  drifts  and  air  currents. 
Both  of  the  beams  exiting  from  the  output  flat  have  the  phase  change 
information  plus  laser  generated  noise.  However  the  phase  change  informa- 
tion is  180°  out  of  phase  between  the  two  beams  whereas  the  laser  noise  is 
in  phase.  Thus  the  output  of  the  differential  amplifier  shows  only  the 
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optical  phase  change  signal.  A lens  is  placed  in  front  of  the  interferometer 
so  that  the  beam  waist  is  in  the  center  of  the  plasma  cell,  increasing  the 
spatial  resolution  of  the  system  to  the  diffraction  limit. 

To  test  the  ability  of  this  system  to  measure  small  phase 
changes,  a CS^  Kerr  cell  was  placed  in  the  interferometer  to  give  a 
measurable  phase  shift.  In  a system  that  used  a Spectra  Physics  166  argon 
ion  laser  with  an  intra-cavity  etalon,  and  an  output  power  of  nearly  1 
watt  at  5145  X,  the  measurements  shown  in  Figure  2.14  were  made.  This 
shows  a 50  KHz  voltage  applied  across  the  Kerr  cell  and  the  resulting 
signal  from  the  differential  amplifier.  This  is  a single  oscilloscope 
trace  recorded  in  real  time  without  the  use  of  any  averaging  techniques. 

It  shows  that  resolutions  of  better  than  0. 1 X are  obtainable  without  any 
trouble.  The  noise  seen  is  due  to  the  laser  beam  shot  noise  and  also  due 
to  very  small  angular  rotations  of  the  laser  beam  due  to  laser  vibrations. 
When  signal  integration  schemes  are  used,  phase  shifts  smaller  than  10  ^ X 
are  readily  observed  with  10  seconds  of  averaging.  Thus  this  interferometer 
is  very  sensitive  and  simple.  To  our  knowledge  this  is  the  most  accurate 
Interferometer  ever  devised  for  the  measurement  of  small  phase  shifts. 

2.5.2  Interferometric  Measurements  on  Plasma  Display  Cells 

A plasma  cell  is  constructed  with  electrodes  of  the  geometry 
shown  in  Figure  2.15.  A visible  laser  beam  can  be  focused  to  a diameter 
of  about  25  pm  over  a 1 cm  path.  The  cell  which  is  made  long  so  as  to 
maximize  the  phase  change  is  moved  relative  to  the  laser  beam  to  determine 
the  phase  changes  at  various  points  along  the  gap. 
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Figure  2.15.  Basic  geometry  of  plasma  cell  electrodes 
used  in  Interferometry  experiments. 
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Figure  2.16  shows  typical  phase  change  results  obtained  with  an 
argon  ion  laser  in  a cell  with  electrodes  with  approximately  the  geometry 
shown  in  Figure  2.15.  The  beam  is  positioned  at  the  center  of  the  discharge 
gap.  Each  discharge  of  the  cell  is  seen  to  initiate  a series  of  strong 
oscillations  that  damp  out  with  time.  These  phase  changes  are  not  typical 
of  electron  density  changes  but  rather  of  neutral  atom  density  changes. 

In  a simple  analysis,  there  are  two  density  changes  that  can  be  measured 
by  measuring  the  phase  change  of  transmitted  light.  The  phase  change  of 
the  electrons,  which  behaves  as  described  by  Equation  (19),  and  the  phase 
change  of  the  neutral  atoms  which  behaves  as  Equation  (18)  where  the 
change  in  the  index,  Ah,  is  proportional  to  the  change  in  density  of  atoms. 
Thus  the  total  phase  change  in  radians  is 

2 ttL ( Cr ( AN  )-l) 

Hot r1 BHLX  (21) 

where  a is  the  proportionality  constant  between  the  atomic  number  density 
Nft,  and  the  index  of  refraction  for  the  particular  gas  atom.  a can  be 
easily  found  from  knowledge  of  the  index  of  refraction  of  the  gas,  n,  at 
a given  temperature  T and  pressure  p by  the  relation  a » ^ where  K is 
Boltzman's  constant.  a is  slightly  dependent  on  wavelength  depending  on 
the  dispersive  characteristics  of  the  gas.  It  is  not  easy  to  separate  the 
two  effects.  The  standard  solution  used  by  most  researchers  is  to  perform 
the  measurement  at  a number  of  wavelengths  and  because  of  the  wavelength 
dependence  of  the  two  terms,  there  is  a possibility  of  separating  them. 
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The  measurement  of  plasma  display  cell  phase  changes  shown  in 
Figure  2.16  shows  rapid  phase  changes  that  are  primarily  due  to  atomic 
density  changes.  The  density  changes  are  due  to  sound  waves  that  are 
initiated  by  the  gas  discharge.  These  sound  waves  travel  back  and  forth 
from  the  anode  to  the  cathode  for  many  cycles  before  they  are  damped  out 
by  absorption.  The  exact  nature  of  these  sound  waves  can  be  determined  by 
noticing  the  phase  changes  for  different  positions  along  the  discharge 
gap.  Figure  2.17  shows  the  cell  discharge  current  and  the  phase  changes 
for  five  differing  positions  along  the  discharge  gap.  The  lowest  trace 
shows  the  discharge  current.  Progressing  upward,  the  next  trace  shows  the 
phase  change  at  the  cathode.  The  third  trace  shows  the  phase  change  54  jm 
from  the  cathode.  The  remaining  traces  are  spaced  by  54  gm,  the  top  trace 
being  at  the  anode.  From  the  calculations  of  Section  2,  it  was  seen  that 
most  of  the  discharge  activity  occurs  near  the  high  field  region  of  the 
cathode.  The  hot  electrons  near  the  cathode  will  of  course  locally  heat 
the  neutral  atoms.  The  more  energetic  neutral  atoms  will  fly  out  of  the 
cathode  region  and  thus  cause  a decrease  in  density.  This  decrease  is 
very  clearly  shown  as  a strong  decrease  in  the  index  of  refraction  at  the 
cathode  at  the  time  of  the  discharge.  The  impulse  of  atoms  flying  out  of 
the  cathode  region  sets  up  a sound  wave  that  is  first  seen  a short  distance 
away  from  the  cathode  as  an  increase  in  index  of  refraction.  The  sound 
wave  travels  through  space  and  time  until  it  hits  the  anode  in  the  top 
trace  and  is  reflected  back  to  the  cathode.  This  bouncing  between  the 
cathode  and  anode  continues  for  many  microseconds  as  shown  in  Figure 
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Figure  2.17. 
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Multiple  exposure  photograph  showing  the  index 
of  refraction  at  5 points  along  the  discharge 
gap,  as  a function  of  time.  The  lowest  trace 
is  the  discharge  current  at  10  ma  per  division. 
The  bottom  index  of  refraction  trace  is  at  the 
cathode  and  the  top  trace  is  at  the  anode. 

The  intermediate  traces  are  spaced  by  54  nm 
intervals.  Horizontal  scale:  .5  p.s  per 

division.  Vertical  scale:  5 K per  division. 


The  question  arises  as  to  how  much  of  the  phase  change  shown  in 
Figures  2.16  and  2.17  is  due  to  electrons.  Clearly  the  rapid  oscillations 
are  due  to  the  neutral  density  changes  since  the  waves  observed  in  Figure 
2.17  travel  at  the  velocity  of  sound.  However,  Figure  2.16  shows  lower 
frequency  components  that  may  be  due  to  electron  density  changes.  To  check 
this  the  interferometer  was  operated  with  a home  built  HeNe  laser  that 
oscillated  simultaneously  on  6328  X and  1.15  ym.  From  Equation  (21),  the 
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ratio  of  the  neutral  change  to  the  electron  change  should  be  ~3  times 
smaller  at  1.15  ym  than  at  6328  X.  Thus  any  electron  phase  changes  that 
are  observable  in  Figure  2.16  should  show  up  much  more  strongly  at  1.15  ym. 

The  sensitivity  of  the  interferometer  with  the  few  milliwatts  coming  from 
the  HeNe  laser,  was  considerably  less  than  that  shown  in  Figure  2.16. 

However  with  a boxcar  integrator,  the  phase  changes  due  to  the  shock 
waves  could  be  resolved  as  good  as  shown  in  Figure  2.16.  As  closely  as 
could  be  measured,  the  phase  changes  in  radians  were  inversely  proportional 
to  wavelength  for  the  two  wavelengths.  Thus  to  the  resolution  of  Figure 
2. 16,  all  of  the  phase  change  is  due  to  neutrals.  This  is  not  too  surprising, 
since  calculations  presented  above  showed  that  an  electron  density  of 
10*  cm  * should  yield  a phase  change  of  only  ~ 10  X.  The  data  shown 
ir.  Figure  2.16  is  scaled  at  ~ 1 X per  division.  Thus  one  could  hardly 
txpect  to  see  an  electron  phase  change  even  at  1.15  ym  with  this  resolution. 
The  neutral  atom  phase  changes  seem  to  amount  to  a 5 X shift.  Thus  there 
could  easily  be  a ratio  of  500  between  the  neutral  and  electron  phase 
shifts.  To  separate  these  two  by  a two  wavelength  experiment  will  require 
an  accuracy  of  the  measurements  much  greater  than  given  by  oscilloscope 
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traces  and  even  the  boxcar  integrators,  due  to  dynamic  range  limitations. 
Measurement  of  electron  densities,  even  with  the  sensitive  interferometer 
setup  described  above,  is  difficult.  This  problem  will  be  discussed 
further  at  the  end  of  this  section. 

2.5.2. 1 Atom  Density  Measurements 

We  first  consider  what  can  be  learned  about  the  discharge  from 

our  ability  to  measure  very  accurately  the  neutral  atom  density  ion  space 

and  time.  A basic  assumption  that  was  used  in  all  calculations  presented 

in  Section  2 was  that  the  neutral  atom  density  was  constant.  This  is  of 

course  an  approximation,  but  how  good  of  an  approximation?  From  Equation 

(21)  and  Figure  2.16,  it  can  be  easily  calculated  that  a 5 X phase  shift 
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represents  only  a 2.0  x 10  cm  change  in  density.  Since  the  atomic 
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density  at  250  torr  and  25  C is  8.84  x 10  cm  , it  is  seen  that  the 
perturbations  are  clearly  very  small,  making  the  constant  density 
approximation  a good  one  for  the  neon  + 0.1%  Xe  mixture.  This  is  probably 
not  true  for  all  gas  mixtures,  however.  An  example  is  shown  for  250  torr 
xenon.  Pure  xenon  is  an  attractive  choice  tor  color  plasma  display  panels 
that  use  the  ultraviolet  light  generated  by  the  xenon  to  excite  phosphors 
[12].  Figure  2.18  shows  the  sound  waves  for  a cell  with  250  torr  xenon. 

The  sound  waves  give  about  a quarter  wavelength  phase  change.  Such  a large 

phase  distorts  the  approximate  linear  relation  between  the  phase  change  and 
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the  output  of  the  interferometer,  however  the  cos  relationship  still 
exists.  This  large  phase  change  is  partially  due  to  the  fact  that  the 
index  of  refraction  at  STP  of  xenon  is  1.000708  [13]  whereas  that  of  neon 
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is  only  1.0000673  [14].  However  the  atomic  density  change  due  to  the 
sound  waves  in  xenon  is  also  very  much  larger  than  that  of  neon.  The  data 
of  Figure  2.18  shows  density  changes  as  great  as  10%  of  the  equilibrium 
density.  The  approximation  of  uniform  neutral  density  starts  to  breakdown 
and  thus  the  effects  of  the  soundwaves  must  be  included  in  discharge 
calculations. 

The  measurements  presented  above  show  only  the  high  frequency 
changes  in  neutral  atom  densities.  There  is,  however,  a steady  state 
change  that  occurs  whether  the  cell  is  on  or  off.  This  density  change  is 
due  to  the  competition  of  the  heating  of  each  discharge  and  the  thermal 
relaxation  or  cooling  of  the  gas.  Thus  this  steady  state  change  is 
established  after  a number  of  discharges  and  is  also  strongly  dependent 
on  the  repetition  rate  of  the  discharges.  To  measure  this,  the  cell  was 
effectively  turned  on  and  off  by  interrupting  the  sustain  generator  for 
many  cycles.  Such  data  is  shown  in  Figure  2.19.  In  this  figure,  the  cell 
which  was  initially  sustained  at  a 50  KHz  rate  was  turned  off  by  stopping 
the  sustainer.  The  neutral  density  quickly  increases  to  its  equilibrium 
value  in  a few  hundred  microseconds.  After  the  equilibrium  value  is  reached, 
there  are  still  small  fluctuations  due  to  soundwaves  traveling  long 
distances  through  the  entire  cell  structure.  After  9 milliseconds  of 
interruption,  the  sustainer  was  again  applied  and  the  discharges  caused  a 
new  steady  state  value  to  be  achieved.  The  wide  band  structure  of  the 
second  half  of  the  trace  is  due  to  all  of  the  soundwaves  bouncing  between 
the  anode  and  cathode  similar  to  Figure  2.16.  Figure  2.20  shows  data  for 
the  same  experiment  but  with  only  a 1 millisecond  sustainer  interruption. 
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Figure  2.19.  Phase  changes  in  Xe  filled  cell  showing  slow  density 
changes  due  to  slow  gas  temperature  changes. 
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The  time  constant  for  the  relaxation  of  the  gas  density  to  equilibrium 
appears  to  be  about  100  |j£ . Figure  2.21  shows  data  for  even  a shorter 
time  scale  where  the  interruption  of  the  sustainer  was  only  about  150  us. 

In  this  case  the  ceil  never  reaches  a steady  state  situation.  The  sound 
waves  reflecting  between  the  cathode  and  anode  are  clearly  visible  here. 

With  each  discharge  the  neutral  atom  density  decreases  slightly  due  to  the 
additional  heating  that  each  discharge  introduces. 

The  thermal  conductivity  of  neon  is  much  greater  than  that  of 
xenon  so  that  the  cell  reaches  steady  state  much  more  quickly.  However 
because  of  this  greater  thermal  conductivity,  less  thermal  energy  can  be 
stored  in  the  gas  and  thus  the  neon  neutral  atom  temperature  much  closer 
to  the  temperature  of  the  glass  wall.  Thus  the  steady  state  density  changes 
in  neon  are  not  as  great  as  those  for  xenon.  These  density  changes  are 
observable,  and  have  characteristics  similar  to  those  of  Figures  2.19,  2.20, 
and  2.21  but  with  much  smaller  phase  changes  and  much  shorter  time  constants. 
However,  due  to  the  low  frequency  drift  characteristics  of  the  interferom- 
eter, it  is  not  possible  to  get  a good  oscilloscope  trace  of  this  data  for 
pure  Ne. 

The  above  mentioned  measurements  indicate  the  nature  of  the 
neutral  density  changes  in  the  gas.  These  density  changes  are  caused 
primarily  by  temperature  changes  in  the  gas.  These  temperature  changes 
in  the  gas  are  indirectly  measurable  by  observing  the  density  changes 
that  they  cause.  For  instance,  in  Figure  2.17,  the  trace  of  the  phase 
change  at  the  cathode  shows  that  at  the  time  of  the  discharge,  the  index 
of  refraction  dropped  sharply.  This  is  a decrease  in  neutral  density 
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Figure  2.21.  Same  conditions  as  in  Figures  2.19  and 
2.20  but  with  time  scale  50  ns/div. 
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directly  attributable  to  a sharp  increase  in  temperature  near  the  cathode 
region.  The  traces  of  Figure  2.17  for  regions  away  from  the  cathode  do 
not  show  this  sharp  drop  at  the  time  of  the  discharge  and  thus  there  is  not 
much  heating.  One  problem  with  this  technique  is  that  the  sound  waves 
also  cause  similar  changes  in  neutral  density.  Thus  it  is  a matter  of 
judgment  and  care  to  come  up  with  the  current  explanation  for  the  observed 
phase  changes.  Nevertheless,  this  is  a useful  qualitative  technique  for 
observing  the  areas  of  heating,  as  will  be  discussed  later  in  this  section. 

A second  problem  with  measuring  the  temperature  change  by  watching 
the  neutral  density  changes  is  that  one  only  finds  out  the  temperature 
changes  and  not  the  actual  temperature.  As  indicated  in  Figures  2.19,  2.20, 
and  2.21,  it  typically  takes  a number  of  discharge  cycles  before  the  gas 
comes  to  an  equilibrium  temperature.  Thus  a technique  is  needed  to  measure 
the  absolute  temperature  of  the  gas.  Fortunately  a simple  technique  exists 
for  doing  this  by  measuring  the  velocity  of  the  sound  waves  bouncing  between 
the  cathode  and  the  anode  as  in  Figure  2.17.  The  velocity  of  sound  in  a 
gas  is  proportional  to  the  square  root  of  temperature  and  very  weakly 
dependent  on  gas  pressure.  The  pressure  dependence  is  too  weak  to  be  of 
concern  for  the  following  calculations.  By  measuring  the  time  for  the 
sound  waves  to  make  many  cathode -anode  transits,  the  velocity  can  be 
accurately  determined  with  knowledge  of  the  cathode -anode  gap  spacing.  To 
make  the  temperature  measurement  less  dependent  on  the  gap  spacing  measure- 
ment, the  sound  wave  velocity  was  measured  at  a number  of  different  sustain 
frequencies.  If  the  sustain  frequency  is  very  low  like  1 KHz,  then  the 
temperature  of  the  gas  will  be  very  close  to  the  ambient  of  the  glass  walls. 
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This  is  true  for  most  all  gases  since  Figure  2.20  shows  that  even  a low 
thermal  conductivity  gas  like  Xe  reaches  an  equilibrium  situation  after  a 
few  hundred  microseconds.  Once  the  velocity  is  measured  at  1 KHz,  the 
sustain  frequency  is  then  increased  to  the  desired  value,  e.g.  50  KHz. 
Because  of  the  dependence  of  velocity  on  temperature,  the  following 
relationship  holds  for  the  two  velocities  and  temperatures: 


T 

T 


H 

L 


2 


(22) 


where  T and  V are  the  temperature  and  velocity  at  the  high  sustain 
H H 

frequency,  and  T and  V are  for  the  low  sustain  frequency.  T is  assumed 

L L Li 

to  be  the  ambient  temperature  of  the  glass  panel.  The  velocities  are  not 
actually  measured,  but  some  convenient  transit  time  which  is  proportional 
to  the  velocity.  This  might  typically  be  the  time  for  10  round  trip 
transists  of  the  sound  wave. 

One  would  expect  that  a cell  filled  with  xenon  should  show  the 
largest  temperature  changes  since  Xe  has  poor  heat  conductivity  compared 
to  the  lighter  inert  gases,  and  also  the  discharge  is  much  more  intense 
in  xenon.  The  transit  times  were  measured  with  a 50  KHz  sustainer  and  a 
1.5  KHz  sustainer.  The  results  showed  that  the  gas  temperature  at  50  KHz 
was  36°C  above  the  temperature  at  1.5  KHz  which  was  assumed  to  be  within 
1°C  of  the  ambient  temperature  of  the  glass.  This  was  with  250  torr  Xe  and 
a 185  ym  gap  spacing.  This  is  the  same  cell  arrangement  used  in  Figures 
2.18,  2.19,  2.20,  and  2.21.  Such  a large  temperature  change  is  consistent 
with  the  ~ 107«  density  changes  shown  in  Figure  2.19  at  a constant  pressure. 


Thus  in  the  case  of  pure  Xe,  the  effects  of  average  gas  temperature  are 
significant  and  must  be  included  in  any  theory  of  operation. 

For  more  standard  gas  mixtures  such  as  Ne  + 0.1%  Xe,  the  temper- 
ature change  is  much  smaller  and  is  not  detectable  with  the  sound  velocity 
technique  because  of  limited  oscilloscope  resolution.  The  resolution 
limits  show  that  the  temperature  increase  cannot  be  greater  than  5°C.  If 
one  assumes  that  the  density  change  such  as  shown  in  Figure  2.19  is 
completely  due  to  a change  in  temperature  at  constant  pressure,  then  the 
temperature  change  in  Ne  + 0.1%  Xe  at  500  torr  and  a gap  spacing  of  185  pm 
is  3.5°C.  Thus  for  typically  used  gas  mixtures  in  plasma  display  panels, 
the  gas  is  always  within  a few  degrees  C of  the  ambient  glass  temperature. 
This  is  not  to  say  that  the  glass  does  not  get  above  the  room  ambient. 

It  is  well-known  that  a fully  lit  plasma  panel  may  rise  10  or  15°C  above 
ambient.  However  the  gas  temperature  will  always  be  within  a few  degrees 
C of  this  warmer  glass  ambient.  Thus  any  model  of  Ne  filled  plasma  panels 
can  justifiably  neglect  the  effects  of  neutral  gas  temperature  change. 

2. 5. 2. 2 Electron  Density  Measurement 

The  laser  interferometry  technique  seems  to  be  an  excellent 

measure  of  the  neutral  atom  density  in  the  plasma  display  panel.  However, 

these  interferometry  experiments  were  performed  in  hopes  of  measuring  the 

electron  density,  a quantity  much  more  important  to  the  operation  of  the 

3 

device.  As  shown  above,  the  neutral  density  changes  may  be  about  10  times 
bigger  than  the  electron  density  changes.  Since  an  electron  density  of 
101  cm  ^ should  give  about  a 10  ^ X phase  change  for  this  experimental 
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setup,  it  is  clear  that  the  problem  is  not  the  signal  to  noise  ratio  of 
the  system  since,  as  shown  in  Chapter  4,  10  ^ X phase  changes  are  above 
the  noise  level  of  this  system  for  reasonable  integration  times.  Thus  in 
theory,  the  electron  density  is  measurable  if  the  electron  and  neutral 
densities  can  be  separated. 

The  standard  technique  for  separating  these  two  densities  is  to 
perform  the  measurement  at  two  different  laser  wavelengths  as  discussed 
above  [15].  From  Equation  (21)  it  is  clearly  desirable  that  the  two 
wavelengths  (\^  and  X^,  X^>\2)  '3e  as  ^ar  aPart  as  possible  in  order  to 
maximize  the  difference  of  the  resulting  phase  changes.  We  can  rewrite 
Equation  (21)  as 


kPtot  * Acpa  + Acpe  <23) 

where  A cca  is  the  phase  change  due  to  neutral  atom  density  changes  and 
A is  the  phase  change  due  to  electron  density  changes.  From  (21)  it  is 
clear  that 

" x (24) 

Acpfi  - EX  (25) 

* 2TTL(a(AN  )-l)  and  E “ -BAn  L.  Since  from  the  above  discussions, 
a e 

2 

x 10  , the  value  of  A can  be  approximated  from  (23)  and  (24)  by 

A ‘ X2AO(tot2 
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where  Acp  is  the  total  phase  change  for  \ 

tot  2 

at  the  two  wavelengths,  one  finds  that 


Xj.  From  (23),  (24),  and  (25) 


E 


X2^, 


totr 


xi*W 


(27) 


Since  Acpg  is  so  small  compared  to  Acpa>  the  question  of  how 
accurate  Acp^^^  must  be  measured  becomes  very  important.  From  (27)  the 
accuracy  of  the  E measurement  is  determined  by  the  accuracy  of  the  numerator 
since  the  wavelengths  in  the  denominator  are  known  to  high  precision.  From 
(23),  (24),  and  (25),  the  numerator  of  (27)  is 


vt; + EV  - Vt2  + EV  (28) 

or  (A  + E\|)  - (A  + EX^  ) . (29) 

2 2 2 2 

Since  A is  very  much  larger  than  E(X^  - \^)  and  E(X^  - Xj)  is  the  sum  of 

the  two  terms  in  the  numerator  of  (27),  a measure  of  the  accuracy  that  the 

values  of  tp  _ must  be  measured  is  found  by  the  ratio 
tot  1 

2 2 
E(Xj  - Xj) 

A 


which  from  (25)  and  (26)  equals 


(30) 
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from  the  wavelength  assignment  convention  — < 1.  Also  from  the  data 

Am  *1 

e i 1 

presented  here  — — < 0.002.  Thus  the  value  of  (30)  is  also  less  than 


0.002.  This  value  represents  the  accuracy  with  which  Acptot  must  be 
known  for  the  error  in  Acf^ot^  to  give  a false  value  of  E equal  to  the  real 
value  of  E in  (27).  Thus  for  greatest  accuracy,  one  should  try  to 
maximize  (30). 


For  signal  to  noise  reasons  an  argon  ion  laser  has  been  chosen 
for  these  experiments.  In  order  to  maximize  (30)  it  is  desirable  to  choose 

and  >.?  as  far  apart  as  possible.  For  the  available  laser  power  the 

Acp 

choice  is  ■ 5145  X and  m 4579  X.  For  — — “ 0.002,  this  gives  a 
-4  ®a 

value  of  4.16  x 10  . Thus  for  the  expected  values  of  Acfe>  and  the 

available  laser,  Aq>tot  must  be  measured  to  an  accuracy  of  better  than  4 
4 

parts  per  10  . Since  from  the  data  of  Figure  2.17,  &Ptot  typically 
5 X,  this  measurement  requires  a resolution  of  at  least  2 x 10  ^ X.  This 
is  clearly  within  the  resolution  of  the  interferometry  system  which  is 
better  than  10  X.  However,  such  resolutions  require  the  use  of  a boxcar 
integrator  for  the  necessary  averaging.  These  instruments  usually  only 

3 

have  a dynamic  range  on  the  order  of  10  . Thus  the  noise  of  the  system 
will  be  determined  by  the  noise  in  the  boxcar  integrator,  even  when  the 
5 X Acptot  input  gives  the  maximum  output  on  the  boxcar  integrator. 

One  approach  to  this  problem  relies  on  the  ability  to  remove  the 
electrons  from  the  gas  volume  by  placing  appropriate  voltage  waveforms  on 
the  cell.  Small  voltage  pulses  placed  on  top  of  the  sustain  waveforms 
have  the  effect  of  sweeping  the  plasma  out  of  the  volume.  The  idea  is  to 
substract  a discharge  cycle  with  plasma  from  a cycle  in  which  the  plasma 
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was  removed.  If  the  only  difference  in  the  two  signals  Is  the  removal  of 
the  plasma,  then  the  plasma  signal  should  show  up  clearly  in  the  difference. 
This  technique  is  now  being  developed. 

2.5.3  Glass  Surface  Properties 

All  of  the  above  experiments  were  directed  at  measuring  properties 
of  the  gas  volume.  Of  equal  importance  are  physical  effects  at  the  glass 
cathode  surface.  A major  question  that  is  still  debated  is  what  proportion 
of  electrons  are  generated  by  photon,  ion  or  metastable  atom  bombardment 
of  the  cathode.  It  has  previously  been  reasoned  that  only  photon 
bombardment  was  important  because  the  speed  of  the  discharge  processes 
was  Taster  than  the  ion  process.  [16,17,18]  This  reasoning  is  false  since 
most  of  the  computer  calculations  performed  in  this  work  used  only  the 
ion  effect  and  achieved  times  as  fast  and  faster  than  the  experimentally 
measured  discharge  times. 

To  resolve  this  issue,  an  experimental  cell  was  made  with  a glass 
cathode.  A pulsed  light  source  was  focused  on  the  glass  cathode  to 
photoemit  a pulse  of  electrons.  With  gas  in  the  cell,  and  the  appropriate 
dc  voltage,  the  electrons  initiate  Townsend  avalanches  which  are  measured 
with  a sensitive  current  detector.  These  avalanches  lead  to  succeeding 
avalanches  by  means  of  the  photon,  ion,  or  metastable  bombardment. 

Because  of  the  grossly  different  speeds  of  these  three  effects,  it  is 
usually  a simple  matter  to  separate  them  in  time.  Figure  2.22  shows  a 
single  trace  of  avalanches  produced  by  a single  light  pulse.  For  exper- 
imental reasons,  the  cell  had  a very  large  gap  spacing  and  thus  a 
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Figure  2.22.  Avalanche  currents  in  a glass  cell  having 
an  electrode  separation  of  2.5  cm  and 
filled  with  5 torr  He  and  0.25  torr  N2* 
Vertical  scale:  5x10“^  Amps  per  division. 

Horizontal  scale:  50  ^s  per  division. 

The  current  increases  downward. 
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necessarily  low  gas  pressure.  The  glass  cathode  was  made  of  common 
window  or  soda  glass.  The  avalanches  were  clearly  regenerated  by  the  ion 
effect  for  this  data.  The  period  of  the  avalanche  group  shown  in  Figure 
2.22  corresponds  to  the  ion  transit  time  across  the  gap.  After  a number 
of  avalanches,  diffusion  processes  cause  the  individual  avalanches  to  be 
obscured.  However  they  continue  until  enough  current  has  flown  to  charge 
the  glass  surface  and  reduce  the  voltage  across  the  discharge  gap.  All 
of  these  experiments  were  performed  with  the  applied  dc  voltage  very 
close  to  the  breakdown  voltage.  No  successive  avalanches  due  to  photon 
bombardment  were  observed.  In  certain  gas  mixtures,  such  as  pure  neon, 
the  ion  processes  were  accompanied  by  a much  slower  growth  process  that 
appeared  to  be  attributable  to  metastable  bombardment. 
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III.  INTERNAL  LOGIC 


3. 1 Introduction 

One  of  the  most  important  implications  of  the  gas  discharge 
studies  is  that  it  should  be  possible  to  enhance  the  internal  logical 
capability  of  the  plasma  display  panel  without  degrading  performance.  This 
section  provides  a background  for  multilevel  internal  decoding,  it  describes 
a new  technique  for  multilevel  writing,  and  it  relates  these  results  to 
complementary  investigations  in  process  at  other  laboratories. 

3.2  Panel  Decoding 

In  a matrix  display  panel  containing  mxn  elements,  a binary 
message  of  (m+n)  bits  must  be  decoded  to  select  a specific  element  for 
writing  or  erasing.  Although  all  of  this  decoding  can  be  performed  rather 
simply  by  circuitry  external  to  the  panel,  this  approach  requires  a control 
circuit  for  each  element,  and  except  for  some  specialized  large  area 
displays  with  modest  information  content  it  is  not  economically  practical. 

A 512  x 512  display,  for  example  would  require  262, 144  control  circuits  and 
at  least  the  same  number  of  connections  to  the  panel. 

If  the  display  elements  have  appropriate  characteristics,  one 
level  of  decoding  can  be  performed  in  the  panel  itself,  and  as  a result, 
the  complexity  of  interconnections  and  external  circuitry  can  be  reduced 
significantly.  In  the  AC  plasma  display  panel,  for  example,  each  display 

element  functions  both  as  a storage  cell  and  as  a logical  "and"  element. 

2 

For  an  n x n panel  the  number  of  addressing  circuits  is  reduced  from  n to 
2n,  and  for  a 512  x 512  panel,  in  particular,  from  262,144  to  1024. 
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The  next  level  of  decoding  is  usually  performed  by  external 
circuitry  which  either  selects  an  active  element  connected  to  the  selected 
row  or  column  electrode,  or  else  provides  a submatrix  of  logical  "and" 
elements  through  which  row  or  column  selection  is  made  by  a smaller  number 
of  active  elements.  For  a 512  x 512  panel  the  "and"  elements  can  be  grouped 
for  each  axis  in  a 32  x 16  submatrix  and  can  be  driven  by  48  active  elements. 

Although  the  circuit  problem  is  manageable,  better  compatibility 
between  panel  and  circuitry  can  be  achieved  if  the  panel  itself  can  provide 
one  more  level  of  decoding.  A 32  x 16  submatrix  of  additional  "and" 
elements  for  each  axis  within  the  panel  requires  only  48  addressing  circuits 
instead  of  512.  If  this  reduction  in  external  circuitry  is  not  offset  by 
increased  panel  complexity,  the  total  cost  of  the  display  panel  and  its 
required  circuitry  can  be  substantially  reduced.  We  describe  several  ways 
of  increasing  the  internal  logical  capability  of  the  display  panel,  with  a 
view  toward  reducing  overall  costs. 

The  essential  structure  which,  with  variations,  is  basic  to  this 
entire  class  of  techniques  replaces  each  row  or  column  electrode  by  a 
pair  of  electrodes,  and  interconnects  these  electrodes  in  two  sets  of 
groups  such  that  selection  of  two  groups  on  an  axis  selects  both  electrodes 
of  only  one  electrode  pair.  For  a 512  x 512  panel  the  electrodes  on  each 
axis  are  divided  into  32  groups  for  one  set,  and  16  groups  in  the  second  or 
orthogonal  set.  The  electrodes  are  driven  by  the  external  circuitry  through 
48  connections.  The  principles  of  this  arrangement  for  a 9 x 9 array  are 
shown  in  Figure  3.1.  In  this  figure  selection  of  group  A and  group  a 
selects  both  electrodes  of  only  one  line  pair. 
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FiRure  3.1.  Electrode  arrangement  for  Four  Input 
Internal  logic. 
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rhe  simplest  approach  to  addressing  in  this  structure  to  provide 
a second  level  of  decoding  was  based  on  the  idea  that  if  a write  signal  were 
applied  simultaneously  to  both  electrodes  in  a pair  the  line  pair  would  be 
fully  selected,  and  at  the  intersection  of  two  fully  selected  line  pairs 
there  would  be  an  initiating  discharge.  If  a write  signal  were  applied  to 
only  one  electrode  of  either  pair  there  would  be  no  initiating  discharge 
Unfortunately  this  configuration  behaves  more  as  a logical  "OR"  circuit 
than  as  a logical  "AND"  circuit  since  the  voltage  on  either  electrode,  or 
both  could  establish  breakdown,  (Because  of  the  wider  electrode  area  provided 
by  two  electrodes,  the  breakdown  would  occur  at  a somewhat  lower  voltage  than 
with  one  electrode,  but  the  margin  for  addressing  would  be  small.) 

A more  promising  variation  of  the  double  electrode  approach  uses 
a metallic  ring  or  a transparent  conducting  dot  to  average  the  voltages 
that  are  applied  to  the  corresponding  electrode  pair.  The  dots  (or  rings) 
are  insulated  from  the  electrodes  by  a thin  dielectric  layer  that  covers 
the  electrodes  and  are  ir Julated  from  the  gas  by  a second  dielectric  layer. 
However,  because  each  element  functions  as  a linear  mixer  there  is  a three 
quarter  select  condition  which  does  not  appear  in  a logical  "and"  circuit 
and  which  could  reduce  operating  margins.  A further  disadvantage  to  the 
system  is  the  more  complicated  structure  that  requires  registration  of  the 
averaging  elements  with  the  intersections  during  assembly.  Nevertheless, 
this  technique  should  function  properly. 

Schernserhorn 's  recently  published  discharge  logic  system  retains 
the  simplicity  of  the  original  structure  and  exploits  the  principle  of 
discharge  coupling  to  discriminate  between  full  selection  and  partial 


67 


I 

I 

I 

I 

i 


I 

| 

i 


selection  [19].  When  two  intersection  line  pairs  are  each  fully  selected 
for  erase,  terminating  discharges  at  each  of  the  four  intersections  end 
the  discharge  sequence  only  at  the  selected  intersections  in  the  element. 

In  the  following  cycles  coupling  from  the  nonterminated  discharges  restores 
the  discharge  sequences  at  the  remaining  intersections - 

Writing  in  this  system,  is  produced  by  erasing,  preceded  and 
followed  by  a state  inversion  which  in  effect  turns  all  "on"  cells  "off" 
and  all  "off"  cells  "on".  In  state  inversion  a bias  voltage  is,  in  effect, 
added  to  the  sustaining  voltage  which  offsets  the  augmenting  effect  of  wall 
voltage.  At  the  same  time  it  increases  the  voltage  across  "off"  cells 
causing  them  to  ignite.  This  system  has  been  demonstrated  by  Scherme rhorn 
at  Owens-Illinois  [19].  For  some  applications  the  aesthetic  problems  of 
int'ersion  and  the  slow  write  cycle  caused  by  the  two  inversion  steps  make 
the  technique  unsuitable. 

3.3  A New  Technique  for  Writing 

The  challenge  in  discharge  logic  is  to  combine  the  simple  erase 
technique  of  Sche rme rhorn 1 s with  an  equally  simple  write  technique,  hope 
for  meeting  this  challenge  is  based  in  part  on  the  fundamental  work 
described  in  section  2 which  indicates  that  ions  and  electrons  formed 
during  a discharge,  or  partially  developed  discharge,  can  remain  in  the 
gap  for  many  microseconds,  that  the  presence  of  these  charges  can  accelerate 
the  development  of  a subsequent  discharge  at  the  same  location,  and  that 
these  charges  can  be  removed  from  the  gap  through  appropriate  pulsed 
voltages.  It  is  also  related  to  work  at  Fujitsu  and  at  the  U of  1 on  shift 
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panels,  and  on  the  discharge  logic  work  at  Owens-Illinois  which  shows  that 
the  presence  of  a discharge  at  one  location  lowers  the  firing  voltage  at 
nearby  locations  [20, 2 1 ] . 

In  one  implementation  of  the  write  technique  pulsed  voltages  are 
applied  to  each  electrode  at  a fully  selected  element.  The  voltages  pulses, 
however,  are  not  applied  simultaneously,  but  instead  are  timed  to  appear 
across  two  of  the  intersecting  electrodes  for  time  and  then  across  the 
remaining  two  intersecting  electrodes  for  time  T^.  A partial  plasma, 
established  during  the  first  interval  , is  drawn  to  the  neighboring 
intersection  during  and  leads  to  a fully  developed  discharge.  The 
discharge  spreads  on  the  next  half  cycle  to  cover  all  four  of  the  electrode 
intersections  at  the  selected  element.  Discharges  should  not  develop  at 
partially  selected  elements. 

In  a variation  of  this  technique  the  initial  select  pulse  height 
is  increased  to  provide  a more  intense  initializing  plasma.  (A  narrow 
pulse  of  high  amplitude  provides  an  intense  discharge  without  developing 
a large  wall  voltage.)  For  a fully  selected  display  element  the  discharge 
spreads  to  the  diagonal  intersection.  For  a partially  selected  element, 
the  tentative  selection  must  be  removed,  either  by  removing  the  plasma  if 
the  wall  voltage  is  small,  or  by  creating  a controlled  discharge  which, 
in  effect,  "erases"  the  portion  of  the  display  element  that  has  been 
tentatively  "written". 

These  techniques  await  the  fabrication  of  experimental  panels 
to  be  fully  investigated.  They  are  discussed  here  because  the  conceptions 
are  results  of  the  fundamental  work,  because  they  appear  to  be  sound,  and 


because  if  successful  they  should  provide  one  method  of  realizing 
significant  cost  savings  through  reduction  of  circuit  complexity. 
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IV.  SPATIAL  GRAY  SCALE 

4.1  Introduction 

In  its  usual  applications  the  plasma  display  panel  displays 
text  and  simple  graphics.  There  are,  however,  many  applications  for  which 
some  degree  of  scale  would  be  desirable.  Probability  density  functions  or 
other  three  dimensional  mathematically  defined  surfaces  could  be  displayed 
using  light  intensity  variations  rather  than  a perspective  drawing. 
Simulations  of  population  growth  or  heat  distribution  could  be  displayed 
in  a similar  manner.  If  a simple  video  input  head  were  developed,  pictures 
could  be  passed  from  one  terminal  location  to  another  in  a facsimile  mode. 

In  fact,  with  a digital  video  input  device  such  as  a CCD,  CID,  or  photodiode 
array,  there  would  be  a direct  geometric  correspondence  between  the  image 
and  display  permitting  direct  digitizing  of  image  point  locations  on  the 
pane  1 . 

In  this  section  we  review  two  methods  of  spatial  gray  scale 
generation,  the  results  of  a hardware  realization  of  one  of  them  (the 
ordered  dither  technique),  and  some  research  on  methods  of  compressing  the 
display  image  data. 

4.2  Dot  Density  Methods 

Unlike  the  CRT  where  each  phosphor  dot  intensity  can  be  varied 
directly,  some  other  form  of  gray  scale  generation  must  be  used  with  the 
bi-level  plasma  panel.  The  current  gray  scale  methods  fall  into  one  of 
two  categories;  pulse  frequency  modulation  or  dot  density  techniques. 
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The  time  division  methods  [22-25],  which  include  the  multi-state, 
the  multi-sustainer , and  duty  cycle  modulation,  vary  the  number  of  light 
output  pulses  per  unit  time  at  each  display  point  to  achieve  gray  scale. 
Because  these  methods  require  extensive  panel  sustainer  modification,  and 
thus  are  not  compatible  with  existing  panel  displays,  our  work  focused  on 
dot  density  techniques. 

Dot  density  gray  scale  relies  on  the  spatial  ratio  of  "on" 
points  to  total  display  points  in  a small  panel  region  to  represent 
intensity  levels.  Since  each  panel  point  is  either  "on"  or  "off"  these 
methods  can  be  displayed  on  any  existing  panel  and  lend  themselves  to 
computer  software  simulation.  These  methods  can  be  classified  into 
adaptive  and  nonadaptive  groups.  In  adaptive  methods  the  decision  to  turn 
a panel  point  "on"  or  "off"  depends  not  only  on  the  intensity  of  the 
original  image  at  that  point  and  possibly  the  position  of  the  point  in  the 
array,  but  also  on  the  character  of  the  image  in  the  neighborhood  of  that 
point.  Nonadaptive  methods  use  only  the  image  intensity  of  the  point  and 
in  some  cases,  the  position  of  the  point.  Because  of  the  complex  nature 
of  adaptive  dot  density  routines,  our  research  has  emphasized  only  non- 
adaptive methods. 

One  of  the  simplest  methods  of  generating  gray  scale  displays 

on  the  panel  uses  a predetermined  pattern  of  1 (ON)  points  in  a NxN  array 

2 

of  panel  points  to  represent  N +1  levels.  The  algorithm  computes  the 
average  of  the  picture  element  (pel)  intensities  in  a NxN  section  of  the 

input  pels  and  maps  the  intensity  onto  the  panel  using  the  nearest  fixed 

2 

pattern  level.  Only  N +1  patterns  are  required  to  generate  any  image,  but 
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since  N pel  intensities  are  averaged,  as  N increases,  the  resolution  of 
the  display  decreases  as  l/N. 

Computer  simulations  of  the  fixed  pattern  method  were  done 
using  two  picture  files.  One  was  a full  gray  scale  picture  with  64 
levels  of  gray  and  the  other  was  a high  contrast  black  on  white  line 
drawing.  The  simulations  show  good  resolution  in  the  5 level  (N“2) 
display,  but  contouring  is  severe  because  of  the  limited  gray  scale. 
Although  sharp  details  are  missing,  17  levels  (N=4)  of  intensity  provide 
a good  compromise  between  gray  scale  and  resolution  for  the  gray  scale 
picture;  but  the  line  drawing  is  greatly  distorted.  The  8x8  pattern  (65 
level)  display  results  in  almost  unrecognizable  reproductions  of  both 
files  due  to  the  large  area  of  video  data  being  averaged.  The  problem  of 
the  fixed  pattern  method  is  the  inability  to  faithfully  reproduce  high 
spatial  frequency  data. 

The  simplest  means  of  generating  a video  display  is  through 
fixed  reference  threshold  detection.  If  the  input  pel  intensity  exceeds 
the  reference,  the  corresponding  panel  point  state  is  set  to  1,  otherwise 
it  is  set  to  0.  For  reproducing  high  contrast  material  such  as  line 
drawings  or  text,  this  technique  works  well,  but  for  other  input  material 
all  of  the  gray  scale  information  is  lost.  However,  if  the  reference  is 
allowed  to  ./ary  randomly  over  the  full  input  intensity  range  for  each  pel 
comparison,  a very  "noisy"  half-tone  representation  results  [26]. 

While  working  with  the  problems  of  TV  video  data  compression, 
Limb  [27]  introduced  a deterministic  signal  which,  when  added  to  the 
digitized  signal  eliminated  some  contouring.  It  was  evident  that  these 
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techniques  could  be  incorporated  on  a bi- level  display  to  generate  gray 
scale.  Known  as  the  ordered  dither  technique  [26],  this  algorithm  compares 
each  scaled  input  pel  intensity  with  a predetermined,  position  dependent 
set  of  threshold  values  to  determine  whether  the  corresponding  display 
point  is  a 1 or  0.  The  size  of  the  so  called  dither  matrix  determines  the 
number  of  distinct  gray  levels  in  the  display  (Figure  4.1).  Because  the  pels 
are  mapped  onto  the  panel  one  to  one,  rapid  intensity  variations  are 
preserved  in  the  display.  Ordered  dither  can  reproduce  regions  of  high 
and  low  spatial  resolution  without  distortion  or  highly  visible  patterns 
[26,29]. 

Computer  simulations  of  ordered  dither  have  shown  this  technique 
to  be  superior  to  the  fixed  pattern  method  in  resolution  and  smooth  gray 
scale  transitions.  Using  ordered  dither  to  generate  the  gray  scale 
photograph  and  the  line  drawing  demonstrates  that  even  the  65  level 
reproductions  display  good  resolution.  It  was  also  found  that  little 
improvement  was  noticed  between  a dither  matrix  of  8x8  (65  levels)  and 
16x16  (256  levels). 

4.3  Hardware  Implementation 

Because  of  its  inherent  simplicity  and  the  encouraging  results, 
ordered  dither  was  selected  for  implementation  in  hardware.  An  interface 
was  designed  to  sample  the  video  output  of  a standard  vidicon  camera  and 
generate  the  necessary  command  words  for  a 16  bit  parallel  input  plasma  panel. 
(Figure  4.2)  The  interface  divides  a TV  video  frame  into  16,  32  bit 
sections.  32  bits  are  sampled  with  a 10  MHz  clock  and  sent  to  the  panel 
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Figure  4 


.1.  The  orderec^dither  matrix.  Each  can 
represent  N +1  intensity  levels. 
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daring  each  horizontal  line  scan  period,  for  a total  display  generation 
time  of  550  ms  [29]. 

Ordered  dither  can  be  generated  in  two  ways.  Either  the  video 
can  be  digitized  and  compared  to  a digital  reference,  or  the  reference 
level  can  be  generated  with  a D/A  converter  and  an  analog  video  comparison 
performed.  The  second  method  was  selected  as  the  least  complex,  and  a 6 
bit  digital  to  analog  converter  was  constructed  using  a resistor  ladder 
network  and  high  speed  transistor  switches  (Figure  4.3).  An  examination 
of  the  required  position  dependent  dither  levels  revealed  that  they  were 
simple  functions  of  X and  Y panel  address.  Logic  circuits  were  designed 
and  built  to  generate  the  proper  binary  codes  for  a 17  and  a 65  level 
ordered  dither  matrix  [29]. 

The  complete  imaging  system  processes  the  input  video  before  it 
i s packed  and  sent  to  the  parallel  panel.  After  viewing  numerous  types 
of  source  images,  the  real  time  system  has  been  shown  to  perform  as  well 
as  the  computer  simulations  predicted  it  should.  The  system  reproduces 
both  full  gray  scale  and  black  on  white  text  well  (Figure  4.4). 

This  video  interface,  constructed  on  five  circuit  cards  using 
standard  TTL  devices,  demonstrates  that  the  hardware  generation  of  video 
displays  on  the  plasma  panel  need  not  be  complex  or  costly. 

4.4  Data  Compression 

The  16  bit  parallel  panel  can  write  or  erase  16  consecutive  Y 
points  each  20  ys,  and  generate  an  entire  512x512  display  in  327  ms, 
which  is  a data  rate  of  800K  bits/second.  Thus,  to  generate  gray  scale 
graphics  in  the  least  possible  time,  transmission  lines  with  bandwidths 
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Figure  4.4. 


An  example  of  an  image  generated 
Digivue  using  the  ordered  dither 


on  a 512-60 
system. 
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much  greater  than  1200  baud  telephone  lines  must  be  used  if  one  is  to 
transmit  each  bit  of  the  display.  However,  since  the  higher  bandwidth 
lines  are  available  only  at  extra  cost,  some  form  of  data  compression  is 
required  to  establish  a compromise  between  line  bandwidth  and  display 
generation  time.  Without  the  ability  to  display  images  in  a reasonably 
short  time  over  low  bandwidth  lines,  and  store  those  images  in  as  little 
memory  as  possible,  few  applications  will  exist  for  panel  graphics. 

Four  display  files  have  been  used  in  these  simulations,  each 
representing  a different  class  of  input.  Two  of  these  files,  LADY  and 
DEMO,  are  full  gray  scale  generated  with  65  level  dither,  one,  IDIOT,  is 
a black  and  white  cartoon  generated  with  dither,  and  the  final  file  is  a 
threshold  detected  version  of  the  cartoon. 

Run  length  coding  has  been  used  with  success  in  transmitting 
threshold  detected  material  rapidly.  [28]  Variable  length  codes  have  been 
shown  to  realize  compressions  of  2 to  4 display  bits  per  transmited  bit, 
but  these  codes  do  not  lend  themselves  to  ordered  dither  because  there  are 
so  many  runs  of  only  a few  bits.  Instead,  a modified  run  length  code  has 
been  devised  known  as  the  4 bit  linear  pattern  code  [29]  that  locates  4 bit 
patterns  in  a line  of  the  display  and  transmits  the  code  word  and  number 
of  consecutive  occurrences.  The  results  using  this  code  are  summarized  in 
Figure  4.5  for  several  input  files.  It  can  be  seen  for  all  of  these  files 
mat  the  compression  ratio  lies  between  1.7  and  2.2.  Thus  this  technique 
would  reduce  transmission  time  from  2.4  minutes  to  a maximum  of  1.4  minutes. 
A maximum  of  32  bits  can  be  represented  by  each  of  the  8 bit  code  words, 
but  even  in  the  worst  case  the  panel  can  keep  up  with  the  data  on  a 1200 
bit/sec  line. 
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PATTERN  SIZE,  N x M 


FILE 

4x4 

4x4^ 

8x4 

2 

8x4 

8x8 

2 

8x8 

4-BLP 

LADY 

2.29 

3.20 

4.00 

5.33 

6.40 

8.00 

2.22 

DEMO 

2. Go 

2.29 

3.56 

4.57 

7.11 

8.00 

1.77 

IDIOT 

2.67 

3.56 

4.57 

6.40 

7.11 

1.81 

IDIOTl 

1.78 

— 

4.78 

— 

6.40 

— 

2.21 

COMPRESSIONS  RATIOS 


1)  Allow  1 bit  error. 

2)  Allow  1,  2,  or  3 bit  errors. 


Figure  4.5.  Compressions  ratios  for  the  4 sample 
files  using  the  N x M pattern  codes 
and  the  4-bit  linear  pattern  code. 
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The  patterns  generated  in  2 dimensions  were  examined  to  determine 
whether  special  characters  could  be  defined  and  displayed  with  character 
codes.  With  such  characters,  the  display  time  would  be  nearly  independent 
of  the  image  composition  once  the  character  set  were  loaded.  The  effective 
compression  ratio  would  be  equal  to  the  character  size  divided  by  the 
character  code  bit  length.  On  a 1200  baud  line,  about  180  6 bit  character 
codes  could  be  sent  per  second  for  a full  screen  display  time  of  82  seconds 
for  4x4  characters  and  21  seconds  for  8x8  characters. 

Investigation  showed  that  very  few  patterns  made  up  most  of  the 
display.  As  Figure  4.6  shows,  the  picture  files  have  10  962  4x4  arrays, 
each  of  which  can  assume  any  one  of  2^=65  536  values,  but  in  two  cases 
fewer  than  75  patterns  exist.  Even  an  8 bit  code  word  reduces  the  storage 
by  50%.  One  can  also  see  how  many  more  patterns  appear  in  the  threshold 
detected  file  because  it  lacks  the  structure  ordered  dither  places  on  each 
array.  Similar  results  were  found  for  8x4  and  8x8  characters. 

Due  to  the  nature  of  dot  density  gray  scale  displays,  a number 
of  isolated  bit  errors  can  be  accepted  in  the  display  before  the  picture 
quality  is  significantly  impared . Therefore,  two  patterns  were  combined 
if  they  differed  by  one  Hamming  distance,  and  the  number  of  patterns  were 
tabulated.  In  all  files  significant  improvement  was  achieved,  reducing 
the  required  number  of  4x4  patterns  to  29  in  one  case.  In  8x4  and  8x8 
characters  2 bit  errors  were  permitted  in  patterns  that  represented  less 
than  .17.  of  the  display,  further  reducing  the  necessary  patterns  to  below 
103  in  all  cases. 

Display  of  the  compressed  picture  files  has  demonstrated  that 
the  errors  can  be  tolerated.  These  results  suggest  that  a set  of  patterns 


PATTERN  SIZE,  N x M 


IDIOTlI  451 


NUMBER  OF  DIFFERENT  PATTERNS 


1)  Allow  1 bit  error. 

2)  Allow  1,  2,  or  3 bit  errors. 


PATTERN  SIZE 


4x4 

8x4 

8x8 

Number  of 
N x M Arrays 
in  the  Files 

10962 

5481 

2709 

TOTAL  NUMBER  OF  ARRAYS 


Figure  4.6.  The  number  of  different  patterns  found  in  each  sample  file 
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and  pattern  codes  could  be  used  to  store  gray  scale  displays  in  perhaps 
207.,  of  the  262K  bits  required  for  a full  panel  bit  by  bit  memory  map. 
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V.  SYSTEMS  RESEARCH 

5.1  Graphic  Display  Terminal/Host  System  Architecture  Studies 

The  design  of  computer  system  architectures  has  been  and  will 
probably  continue  to  be  closely  governed  by  the  nature  of  the  tasks  that 
the  system  is  required  to  perform.  It  seems  clear  that  no  single  optimum 
arrangement  of  components  and  subsystems  will  provide  for  the  wide  variety 
of  needs  of  the  Navy.  Although  there  will  remain  differences  in  overall 
architecture  there  are  some  subsystems  which  are  required  to  exhibit  similar 
characteristics  from  system  to  system.  Important  examples  are  the  terminal 
and  communication  subsystems  which  will  be  used  in  future  ground  based 
computer-based  information  systems  and  computer-based  instrumentation  and 
monitoring  systems  in  aircraft,  submarines,  and  ships. 

The  primary  purpose  of  a computer-based  information  system  is  to 
provide  a user  with  rapid  access  to  stored  or  computer  generated  information 
usually  available  from  some  central  source,  algorithm,  or  data  bank.  There 
are  a considerable  number  of  large  information  systems  emerging  at  the 
present  time  through  the  Navy.  These  include  ARPANET,  Lockheed  Information 
Services,  systems  which  will  supply  maintenance  and  calibration  support, 
both  ground  based  and  on-board  education  and  consulting,  as  well  as  real 
time  instrumentation  and  monitoring.  In  each  of  these  systems  however  there 
is  a common  configuration  emerging,  forced  both  by  technology  and  economics. 
It  is  basically  that  of  a host  system  communicating  in  a conversational  mode 
to  a rather  sophisticated  intelligent  graphics  terminal  over  voice  grade 
communication  lines.  In  the  large  number  of  general  systems  one  does  not 
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expect  this  basic  configuration  to  change  radically  even  though  the  configu- 
ration of  the  host  and  the  terminal  may  vary  considerably  because  of  changes 
in  display  technology,  low  cost  mass  storage,  mini  and  microcomputers  etc. 

It  is  not  expected  that  the  cost  of  bandwidth  between  host  and  terminal  will 
drop  in  the  near  future,  in  fact  the  opposite  will  probably  be  the  case. 

As  a result,  we  have  felt  that  it  is  important  to  continue  to  investigate 
ways  in  which  the  "apparent"  bandwidth  between  the  user  and  the  host  can  be 
optimized.  By  this  we  mean  ways  in  which  the  human  interaction  with  the 
information  can  take  place  in  a natural  manner.  We  would  like  to  emphasize 
that  what  is  important  in  any  information  system  is  the  efficient  and 
economic  transfer  of  information  to  and  from  the  human  requesting  service. 

It  is  important  that  system  designers  and  engineers  recognize  that,  in  the 
end,  it  is  overall  efficiency  of  information  transfer  that  is  the  important 
parameter  that  needs  to  direct  the  system  architecture.  For  example  the 
increased  productivity  of  users  who  can,  in  real  time,  manipulate  models, 
monitor  activity,  access  data  and  consult  with  an  information  bank  through 
the  use  of  an  interactive  graphics  terminal  is  quite  remarkable. 

One  basic  result  of  our  work  in  this  area  has  been  the  development 
of  a systems  model  which  takes  the  form  shown  in  Figure  5.1  and  Figure  5.2. 
Figure  5.1  illustrates  the  arrangement  which  we  feel  will  emerge  in  the  case 
of  large  ground  based  systems  even  though  remote  terminals  may  be  aboard 
ship.  Figure  5.2  illustrates  the  general  arrangement  that  seems  to  be  emerging 
in  the  case  of  local  instrumentation  systems  such  as  those  used  in  future 
aircraft,  submarines,  and  ships.  In  each  case,  it  seems  clear  that 
significant  cost  performance  benefits  can  be  achieved  by  using  terminal 
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Figure  5.1.  Architecture  of  a large  computer-based  information  system. 
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Control  Data  Corp.  Cyber  70  computer  that  services  more  than  500  interactive 
graphics  terminals  over  voice  grade  telephone  lines.  In  addition  to  the 
PLATO  IV  terminal  simulator,  other  terminal  software  systems  are  being 
developed,  including  an  ASCII  Teletype  and  an  IBM  3270.  Throughout  the 
design  of  both  the  hardware  and  software  of  this  terminal,  an  effort  has 
been  made  to  retain  maximum  flexibility  while  observing  the  cost -performance 
considerations  required  in  the  real  user  environment.  Following  is  a 
detailed  summary  of  the  hardware  and  software  architecture  of  the  intelligen 
terminal.  Also  included  is  a description  of  some  methods  of  improving  the 
performance  of  a large  scale  timesharing  system  which  have  been  developed 
using  the  intelligent  terminal  on  the  PLATO  IV  system. 

5.2.1  Hardware  Architecture 

Figure  5.3  shows  the  hardware  arrangement  used  in  the  experimental 
intelligent  terminal.  The  PDP  11/10  mini-computer  is  an  unmodified  8K 
core  configuration  with  four  "small  peripheral  controller  slots"  available. 
The  small  peripheral  controller  slots  were  used  to  house  the  various 
terminal  subsystem  interfaces,  a ROM  bootstrap,  the  display  interface, 
the  communications  interface,  and  the  external  equipment  interface. 

The  display  used  in  the  terminal  is  a modified  Digivue  512-60 
display  memory  unit.  Although  the  display  format  is  similar  to  that  used 
in  the  standard  PLATO  terminal,  this  unit  has  been  upgraded  to  operate  at 
character  writing  rates  of  up  to  6000  characters  per  second  by  means  of  a 
parallel  writing  technique  [31].  This  feature  is  realized  by  modifying 
the  display  address  circuitry. 


The  present  communications  interface  accommodates  PLATO 
input/out*-  it  format  (20  bit  words  in;  10  bit  words  out)  and  ASCII/EIARS232 
signals  via  program  control.  Work  is  presently  in  progress  to  allow  the 
format  of  the  communications  interface  signals  to  be  modified  under  program 
control  in  order  to  provide  automatic  access  to  systems  which  require 
different  protocols.  This  approach,  while  it  is  somewhat  similar  to  the 
DEC  DL11E  communications  interface,  has  an  expanded  maximum  word  size  to 
provide  for  added  format  variability. 

The  external  equipment  interface,  the  DRll-C,  supports  user  input 
devices  such  as  keysets  of  various  types,  touch  input  systems,  and  in 
later  versions,  manual  manipulation  techniques  such  as  light  pens,  joy 
sticks,  etc.  The  external  equipment  interface  is  also  available  for  use 
as  an  input/output  port  for  user  controlled  devices  such  as  laboratory 
equipment,  educational  aids,  and  office  equipment. 

As  the  library  of  terminal  software  grows  and  as  additional 
features  of  the  intelligent  terminal  (such  as  expanded  character  sets  and 
voice  input)  are  implemented,  auxiliary  forms  of  low  cost  storage  will  be 
needed.  For  example,  to  facilitate  a rather  sophisticated  bootstrap 
procedure,  a 256,  16  bit  word  ROM  has  been  added  to  the  terminal.  It  is 
expected  that  as  character  sets,  conversion  tables,  and  various  loaders 
are  developed,  this  memory  will  be  expanded.  For  the  larger  amounts  of 
auxiliary  storage  that  will  be  needed,  configurations  of  solid  state  shift 
register  memory  and  floppy  disk  are  being  evaluated.  For  example,  an  8K 
MOS  shift  register  memory  is  being  designed  and  tested  that  exhibits  a 
worse  case  access  time  of  1 ms  and  a transfer  rate  of  1 MHz.  An  important 
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aspect  of  this  research  is  the  development  of  optimum  memory  hierarchies 
which  will  improve  the  cost -performance  characteristics  of  an  intelligent 
terminal. 

5.2.2  Terminal  Software 

To  best  accommodate  different  host  systems,  the  terminal  software 
is  arranged  as  a single  resident  terminal  monitor,  with  a number  of  host 
dependent,  and  user  dependent  tasks  as  shown  in  Figure  5.4.  The  combination 
of  the  monitor  and  a host  dependent  task  comprise  the  nucleus  of  the 
terminal  simulator  for  each  particular  host.  All  other  routines  are  defined 
as  user  dependent  tasks.  Currently,  the  intelligent  terminal  is  capable 
of  two  modes  of  operation:  as  a PLATO  IV  student  terminal,  and  as  an  ASCII 

teletype.  Other  configurations  under  consideration  are:  a graphics 

terminal  compatible  with  the  ARPANET  graphics  protocol,  an  IBM  3272  for 
access  to  information  and  retrieval  systems,  and  an  extended  stand-alone 
mode  such  as  single  user  BASIC  with  added  graphics  features. 

The  terminal  monitor  contains  all  the  software  drivers  for  the 
plasma  panel.  Including  the  character  generator,  a local  mode  keyset 
handler  with  some  single  line  editing  capabilities,  the  file  handling 
system,  and  all  user  information.  It  is  capable  of  a limited  stand-alone 
mode  of  operation  and  serves  as  a transfer  point  from  one  mode  of  operation 
to  another.  A user  may  "escape"  to  the  monitor  at  any  time  while  using  the 
intelligent  terminal.  A new  host  dependent  task  can  then  be  called,  and 
the  intelligent  terminal  will  automatically  switch  to  the  new  mode, 
establishing  communication  with  a different  host  system  If  necessary.  Most 
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terminal  monitor  subroutines  may  also  be  used  by  either  host  dependent  or 
user  dependent  tasks. 

Each  host  dependent  task  must  interpret  data  from  the  communica- 
tions interface  such  that  the  information  is  either  used  to  call  subroutines 
resident  in  the  terminal  monitor,  passed  to  user  dependent  subroutines,  or 
used  within  the  host  dependent  task.  It  may  handle  all  the  user  input 
devices  directly,  or  through  the  monitor  subroutines.  For  most  applications, 
only  one  host  dependent  task  need  be  resident.  However,  if  inter-host 
communication  is  desired,  the  tasks  for  each  host  must  be  available.  It 
should  be  noted  that  while  the  possibility  of  attaching  several  displays 
and  input  devices  so  as  to  allow  several  users  to  work  simultaneously 
through  one  mini-processor  exists,  this  problem  has  not  been  studied  in 
this  project  to  date. 

The  most  developed  host  dependent  task  is  the  PLATO  IV  terminal 
simulator.  This  task  contains  all  the  features  described  by  Stifle  [30], 
plus  several  new  features  such  as  fast  area  erasure  or  "block  erase",  and 
the  facilities  for  communication  between  TUTOR  [32]  programs  and  terminal 
resident  assembly  language  subroutines. 

Another  mode  of  operation  which  is  now  available  is  the  simulation 
of  an  ASCII  teletype  with  display.  This  has  been  used  to  communicate  with 
a Burroughs  6700,  and  a DEC  System  10.  Inter-host  communication  has  been 
achieved  by  formating  information  on  the  P1AT0  IV  system  and  sending  it  to 
the  B6700  for  extended  processing.  Plans  have  been  made  to  communicate 
with  other  large  scale  information  and  retrieval  systems  through  the 
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intelligent  terminal,  and  to  send  the  selected  information  to  the  PLATO  IV 
system  to  be  displayed  as  part  of  information  presented  as  lesson  material. 

User  dependent  tasks  are  those  which  provide  additional  features 
to  the  basic  terminal  structure.  These  include  terminal  resident  sub- 
routines which  communicate  with  programs  executing  on  the  host  system, 
such  as  those  which  enhance  graphics,  and  software  drivers  that  control 
special  peripheral  equipment.  The  mini-computer  can  be  configured  as  a 
process  control  machine,  yet  simultaneously  retain  the  capability  of 
accepting  input  and  transmitting  data  as  a graphics  terminal.  Since  most 
of  the  initial  work  on  this  terminal  has  been  on  the  PLATO  IV  system, 
nearly  all  the  user  dependent  tasks  now  working  address  themselves  to  the 
problems  of  a large  scale  timesharing  system  with  a relatively  low  bandwidth 
between  the  central  processor  and  the  terminals.  Some  of  these  tasks  are 
described  later  in  this  paper. 

5.3  Specific  Studies  Relating  to  the  PLATO  IV  System 

Throughout  the  development  of  the  experimental  terminal  system, 
studying  the  relationship  between  the  host  computer  and  the  intelligent 
terminal  processor  has  been  an  important  aspect  of  the  research.  For 
example,  rather  than  develop  a local  operating  system  on  the  PDP  11/10, 
the  facilities  of  the  host  system  were  used  whenever  possible.  A cross 
assembler  for  the  PDP11  assembler  language  (PAL)  has  been  written  in  TUTOR 
[32]  the  educational  language  for  PLATO.  This  cross  assembler  runs  as  a 
normal  TUTOR  lesson.  The  binary  output  of  this  assembler  is  saved  on  disk 
in  the  host  system.  These  binaries  can  then  be  inspected,  linked,  and 


sent  to  the  terminal  simulator  from  another  TUTOR  lesson.  A protocol  has 
been  established  so  the  user  may  "log-in"  any  binary  which  is  loaded  in 
this  manner,  identifying  it  by  a six  character  name.  The  resident  program 
can  then  be  called  by  name  from  a TUTOR  program.  Source  code  for  assembly 
is  stored  in  the  same  manner  as  TUTOR  source  code,  so  that  all  of  the  PLATO 
editing  and  filing  systems  are  available  for  use. 

Another  consideration  was  the  bootstrapping  procedure.  A method 
was  desired  that  did  not  require  any  additional  peripherals;  so  again  the 
host  system  facilities  have  been  used.  A ROM-held  bootstrap  program  sends 
a signal  to  PLATO  indicating  that  the  terminal  is  ready  to  be  initialized. 
P1ATO  then  responds  by  sending  the  absolute  loader  in  bootstrap  format, 
and  then  the  monitor  and  PIATO  IV  terminal  simulator  in  absolute  loader 
format.  Control  is  apssed  from  the  bootstrap  to  the  loader  to  the  monitor 
without  user  intervention.  Once  the  monitor  is  running,  the  user  can 
request  which  type  of  terminal  is  desired.  The  appropriate  host  dependent 
task  is  loaded  and  control  is  automatically  transferred  to  that  task. 
Currently,  all  files  are  kept  on  the  PLATO  system.  However  similar 
methods  of  bootstrapping  from  other  host  systems  or  from  local  storage 
could  be  devised. 

Another  basic  problem  of  using  a system  based  on  a large, 
centrally  located  computer  is  the  bandvldth  limitations  of  the  voice  grade 
conmunlcation  lines  which  form  the  connection  between  the  processing  power 
and  the  graphic  displays.  A relatively  small  amount  of  local  processing 
capability  can  effectively  Increase  the  bandwidth  in  several  ways.  For 
example,  it  is  a well-known  fact  that  the  frequency  of  occurrence  of 
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symbols  and  words  In  s language  can  be  determined,  and  that  this  relation- 
ship is  mostly  Independent  of  context.  This  relation  can  be  used  to  enhance 
the  performance  of  a terminal  even  in  a fixed  bandwidth  system  by  using 
more  efficient  coding  schemes.  The  use  of  more  efficient  codes  has  been 
generally  avoided  in  terminals  because  of  the  increased  cost  and  complexity. 
With  the  development  of  low-cost  micro-  and  mini -processors  and  memory, 
use  of  more  efficient  codes  seems  possible  with  little  Increase  in  coat. 

An  investigation  of  this  concept  has  been  initiated  using  the  Intelligent 
terminal  system. 

Another  use  of  a resident  processor  is  local  curve  generation. 
Circles  and  other  conic  sections  have  been  produced  with  significant 
Improvement  in  generation  speed.  In  addition,  the  extensive  use  of  the 
programmable  characters  to  create  multicharacter  pictures  has  prompted 
the  following  work.  In  any  picture  of  more  than  two  characters,  or  roughly 
3/4"  x 3/4",  the  individual  characters  are  plotted  so  slowly  as  to  be 
noticeable.  This  problem  can  be  alleviated  by  sending  encoded  information 
to  the  intelligent  terminal  and  allowing  the  mini-processor  to  plot  the 
entire  figure  at  the  speed  of  the  parallel  panel.  A large,  block  alphabetic 
character  set,  each  character  of  which  is  made  up  of  two  to  six  programmable 
characters,  can  now  be  written  at  the  same  speed  as  the  8x16  dot  characters. 
This  is  an  increase  of  greater  than  six  times  normally  available  speeds, 
sinca  extra  positioning  information  must  normally  be  sent  with  the 
characters. 
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5.4  Conclusions  of  System  Studies 

The  recent  advances  in  micro-processor,  memory,  and  display 
technology  suggest  that  intelligent  terminal  architectures,  which  exhibit 
interactive,  graphic  capability  and  multi-host  compatibility,  will  become 
economically  viable  to  users  of  the  emerging  computer-based  Information 
system  services.  The  work  described  in  this  paper  has  dealt  with  the 
development  of  an  experimental  facility  which  can  be  used  to  model  and 
test  both  intelligent  terminal  architecture  concepts  and  terminal  hardware 
techniques.  This  facility  has  already  been  used  to  develop  and  demonstrate 
numerous  terminal  improvements  and  host/terminal  software  concepts.  These 
include  more  efficient  use  of  the  communication  bandwidth  between  the 
terminal  and  the  host;  automatic  or  programmed  controlled  access  to  a 
variety  of  host  systems  which  require  different  communication  protocol; 
and  Improved  task  performance  as  a result  of  coordination  between  a 
centrally  located  large  scale  host  system  and  a local  mini-computer-based 
terminal  system. 
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